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^ Abstract. This paper presents the results of a high resolution radio imaging survey of all known (96) low- 

, luminosity active galactic nuclei (LLAGNs) at D < 19 Mpc. We first report new 2 cm (150 mas resolution using 

f — ' the VLA) and 6 cm (2 mas resolution using the VLBA) radio observations of the previously unobserved nuclei 

I in our samples and then present results on the complete survey. We find that almost half of all LINERs and low- 

■ luminosity Seyferts have fiat-spectrum radio cores when observed at 150 mas resolution. Higher (2 mas) resolution 

' observations of a fiux-limited subsample have provided a 100% (16 of 16) detection rate of pc-scale radio cores, 

, with implied brightness temperatures ^ 10^ K. The five LLAGNs with the highest core radio fiuxes also have 

^ pc-scale 'jets.' Compact radio cores are almost exclusively found in massive ellipticals and in type 1 nuclei (i.e. 

nuclei with broad Ha emission). Only a few 'transition' nuclei have compact radio cores; those detected in the 
radio have optical emission-line diagnostic ratios close to those of LINERs/Seyferts. This indicates that some 
transition nuclei are truly composite Seyfert/LINER-|-HII region nuclei, with the radio core power depending on 
^ I the strength of the former component. The core radio power is correlated with the nuclear optical 'broad' Ha 

• • ■ luminosity, the nuclear optical 'narrow' emission-line luminosity and width, and with the galaxy luminosity. In 

these correlations LLAGNs fall close to the low-luminosity extrapolations of more powerful AGNs. The scalings 
suggest that many of the radio-non-detected LLAGNs are simply lower power versions of the radio-detected 
^ > ^ ■ LLAGNs. The ratio of core radio power to nuclear optical emission-line luminosity increases with increasing bulge 

' luminosity for all LLAGNs. Also, there is evidence that the luminosity of the disk component of the galaxy is 

correlated with the nuclear emission- line luminosity (but not the core radio power). About half of all LLAGNs 
with multiple epoch data show significant inter-year radio variability. Investigation of a sample of ~ 150 nearby 
bright galaxies, most of them LLAGNs, shows that the nuclear (<150 mas size) radio power is strongly correlated 
with both the black hole mass and the galaxy bulge luminosity; linear regression fits to all ~ 150 galaxies give: 
log P2cm = 1.31(±0.16) log Mmdo -1-8.77 and log P2cm = 1.89(±0.21) log LB(bulge) -0.17. Low accretion rates 
(< 10~^ — 10~^ of the Eddington rate) are implied in both advection- and jet-type models. In brief, all evidence 
points towards the presence of accreting massive black holes in a large fraction, perhaps all, of LLAGNs, with the 
nuclear radio emission originating in either the accretion inflow onto the massive black hole or from jets launched 
by this black hole - accretion disk system. 

Key words, accretion, accretion disks — galaxies: active — galaxies: jets — galaxies: nuclei — radio continuum: 
galaxies — surveys 
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1. Introduction 

The debate on the power source of low-luminosity ac- 
tive galactic nuclei (LLAGNs, i.e. low-luminosity Seyferts, 
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LINERs, and "transition" nuclei [nuclei with spectra in- 
termediate between those of LINERs and HII regions]) 
is a continuing one. Their low emission-line luminosities 
(Lho < 10^° erg s~^ by definition; Ho ct al., 1997a , here- 
after H97a) can be modeled in terms of photoionization by 
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hot, ycung stars (Terlevich & Melnick, 1985; Filippenko 
& Tcrl ;vich, 1992; Shields, 1992|), by coUisional ionization 



fc Wilson, 1989|; [Vila et al., 1990| ; [Sadler et al., 1995| ) 



though Ulvestad & Ho (2001a) have pointed out that 



in shocks (Koski & Osterbrock, 1976; Fosbury et al., 1978; 
Heckman, 1980; Dopita fc Sutherland, 1995| ) or by aging 



starbursts ( Alonso-Herrero et al., 1999| ). 

On the other hand, evidence has been accumulating 
that at least some fraction of LLAGNs share characteris- 
tics in common with their more powerful counterparts - 
radio galaxies and powerful Seyfert galaxies. These simi- 
larities include the presence of compact nuclear radio cores 



(Heckman, 1980), water vapor megamascrs (Braatz ct al 
1997), nuclear point-like UV sources ( Maoz et al., 1995|; 



Barth et al., 1998p , broad Ha lines ( po et al., 1997c| , here 
after H97c), and broader Ha lines in polarized emission 
than in total emission ( Barth et al., 1999 ). If LLAGNs are 
truly scaled down AGNs then the challenge is to explain 
their much lower accretion luminosities. This requires ei- 
ther very low accretion rates 10~® of the Eddington 
accretion rate, LEdd) or radiative efficiencies (the ratio of 
radiated energy to accreted mass) much lower than the 
typical value o f ^ 10% (e.g. Chapter 7.8 of Frank, King, 
& Rail e, 1995) assumed for powerful AGNs. 

One well-known property of some powerful AGNs is 
a compact (sub-parsec), flat-spectrum (usually defined as 
a > —0.5, Sj/ oc v°') nuclear radio source, usually in- 
terpreted as the synchrotron self-absorbed base of the jet 
which fuels larger-scale radio emission. Astrophysical jets 
are known to be produced in systems undergoing accre- 



low-luminosity Seyferts have a higher incidence of flat- 
spectrum cores than more powerful Seyferts. 

How does one distinguish accretion-powered LLAGNs 
from LLAGNs powered by hot stars or supernova shocks? 
Broad Ha lines and bright unresolved optical or UV 
sources are ambiguous indicators because they can all be 
produced in starburst models (see Terlevich et al., 1992 ). 
Searches for broader polarized Ha emission are currently 
feasible in only a few of the brightest LLAGNs (Barth et 
1999). X-ray emission from the nucleus can be con- 



al 



tion onto a compact object (e.g. Pringle, 1993; Blandford 
1993|)~^ such compact radio sources in galactic nuclei may 
reasonably be considered a signature of an AGN. Much 



theoretical work (e.g. 


Begelman, Blandford, & Rees, 1984|; 


Lovelace & Romanova, 1996; 


Falcke & Biermann, 1999) 



has been devoted to this disk-jet relationship in the case of 
galactic nuclei and it has been suggested that scaled-down 
versions of A GN jets can produce flat-sp ectrum radio cores 
in LLAGNs ( Falcke fc Biermann, 1999| ). Compact nuclear 



radio emission with a flat to inverted spectrum is also ex- 
pected fr om the accretion inflo w in advection-dominated 
(ADAF; Narayan et al., 1998 ) or convection-dominated 
(CDAF; Narayan et al., 200C) accretion flows, possible 



f orms of accretion onto a black hole at low accretion rates 
(Rees al., lat-spectrum radio sources can also 



fused with sources in the bulge or disk of the galaxy, even 
at Chandra^s sub arcsecond resolution. Further, all these 
indicators may be affected by viewing geometry, obscura- 
tion, and the signal-to-noise of the observations, the last 
problem being exacerbated by the low optical, UV, and 
X-ray luminosities of LLAGNs and the need to subtract 
the starlight. The radio regime, however, offers several ad- 
vantages. Gigahertz (cm wave) radiation does not suffer 
the obscuration that affects the UV to infra-red. Also, at 
tens of gigahertz the problems of free-free absorption can 
be avoided in most cases. Finally, high resolution, high 
sensitivity radio maps can be routinely made with an in- 
vestment of les s than an hour per sour ce at the Very Large 
ArrayQ (VLA; Thompson et al., 198C) and the Very Long 
Baseline Array (VLBA; [Napier et al., 1994| ) . At their res- 
olutions of ~ 100 milli-arcsec (mas) and ~ 1 mas, respec- 
tively, it is easy to pick out the AGN, since any other radio 
emission from the galaxy is usually resolved out. 

Closely related to these theoretical and observational 
issues are the increasing number of accurate mass deter- 
minations for "massive dark objects" (MDOs) in nearby 
galactic nuclei. Indeed, dynamic signatures of nuclear 
black holes are being found in almost every nearby galaxy 
studied. Mass estimates from high-resolution stellar-, 
maser-, or gas-dynamics are probably accurate to a fac- 
tor of two ( [Richstone et al., 1998| ; [Gebhardt et al., 200Q[ ). 
Slightly less reliable (factor ^ 2-5) estimates of Mmdo for 
large samples are possible using the tight correlation be- 
tween central stellar velocity dispersion of the galaxy bulge 
(cTr) and Mmdo ([Merritt fc Ferrarese, 200l[; [Gebhardt 



result through thermal emission trom ionized gas m nor- 
mal H II regions or through free-free absorption of non- 
thermal radio emission, a process which probably occurs 
in compact nuclear starbursts ( Condon ct al., 1991 ). The 
brightness temperature, Tb, in such starbursts is limited 
to log [Tb (K)] ^ 5 ( [Condon et al., 1991] ). Thus it is nec- 
essary to show that Tb exceeds this limit before accretion 
onto a black hole can be claimed as the power source. 



From the observational perspective, Heckman (1980) 
showed that LINER nuclei tend to be associated with a 
compact radio source, and compact, flat-spectrum radio 
cores are known to be present in many 'normal' E/SO 



et al., 2000). In this paper we assume that the few x 
10^5-9 Mq MDOs traced by dynamic methods are black 
holes, and interchangeably use the terms MDO and black 
hole. With the new MDO results the degeneracy between 
black hole mass and accretion rate when accounting for 
radiated luminosity is removed. Given AGN spectral en- 
ergy densities and black hole masses in large samples of 
nearby AGNs, one can attempt to isolate factors which 
determine accretion rates and luminosities. 

A weak correlation between the nuclear centimeter ra- 
dio luminosity and the mass of the nuclear black hole 



galaxies (Sadler et al., 1989; Wrobel fc Heeschen, 1991 



Slee et al., 1994). Flat-spectrum radio cores are, however. 



uncommon in normal spirals or Seyfert galaxies ( Ulvestad 



^ The VLA and VLBA are operated by the National 
Radio Astronomy Observatory, a facility of the National 
Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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has been found by several authors, e.g. Franceschini et al 
(1998|)] |Yi fc Boughn (1999| ) and |Di Matteo et al. (200l| ) 
This result has been used as evidence for the existence of 
an ADAF-like inflow since, in such a flow, the radio lumi- 
nosity is predicted to scale with the black hole mass and 
the accretion rate (e.g. Mahadevan, 1997). Two signifi- 
cant weaknesses here are the low resolutions of the radio 
data, which allows significant contamination from other 
processes in the bulge, and the small sample sizes. Our 
150 mas resolution radio observations of a large number 
of nearby LLAGNs considerably increase the number of 
LLAGNs with reliable black hole mass estimates and high 
resolution radio observations, allowing a better test of the 
relationship between the two quantities. 

In the following sections we first define the samples 
used in the paper (Sect. ||), and report on new VLA 
and VLBA observations which complete our survey of 
the distance-limited (D < 19 Mpc) sample of LLAGNs 
(Sects. ||and^). Observational results on radio- variability 
in LLAGNs are also presented in Sect. ^. In Sects. || to || 
we combine the new observations with those reported in 
the earlier papers in this series, and present comprehensive 
results on the radio properties of LLAGNs in the distance- 
limited (D < 19 Mpc) sample, and their implications for 
the central engines in LLAGNs. In Sect. ^ we study the 
relationship between nuclear radio power and black hole 
mass, both for the distance-limited LLAGN sample and 
for the larger sample of nearby bright galaxies (most of 
them LLAGNs) defined in Sect. The results are fol- 
lowed by a brief discussion (Sect. |l^) on the nature of the 
central engines in LLAGNs. 

2. Sample Selection 

2.1. The LLAGN Sample 

The results in this paper, and the new observations re- 
ported here, are based on LLAGNs selected from the 
Palomar spectroscopic survey of all (^ 470) northern 
galaxies with Bt < 12.5 mag ( Ho et al., 1995| ). Almost 
half of all galaxies in the Palomar survey host some form 
of LLAGN (Ho et al., 1997b). Spectroscopic parameters 
(including activity classification) of 403 galaxies in the 
Palomar sample have been presented in H97a, and from 
these we have chosen to study the 96 nearest (D < 19 Mpc) 
LLAGNs - the "distance-hmited LLAGN sample" - at 
2 cm with the VLA at ~ 0'.'15 resolution. The distribution 
of galaxy distance for the distance-limited LLAGN sam- 
ple is shown in Fig. Clearly, biases arising from the 
distance distribution are expected to be small within this 
sample. Observations at 2 cm of 33 LLAGNs in the sam- 
ple (take n in 1996 October) were presented in Nagar et| 
al. (2000|) . New 2 cm observations of 72 LLAGNs in the 



J VLA 



sample are presented here (10 LLAGNs are in common 
between this work and Nagar et al. 2000, and 1 LLAGN 
[IC 356] was not observed). 

Our VLA observations of the 95 LLAGNs resulted in 
the detection of 31 nuclei. Sixteen of these have core flux 



> 2.7 mJy, and all sixteen cores have a flat radio 
spectrum^ (a > —0.3, S^, (x j/"). These sixteen objects 
- the "distance- and radio-flux- limited LLAGN sample", 
listed in Table |l[ were chosen for further study at multi- 
ple frequencies with the VLA and at higher resolution with 
the VLBA. Multifrequency VLA observations of the nu- 
clei in Table (Nagar, Wilson, & Falcke 2001) and VLBA 
observations of roughly half of the nuclei in Table |^ have 
already been published (see column 5), while VLBA ob- 
servations of the remaining seven nuclei in Table |^ are 
presented in the present paper. 



2.2. The "MDO" Sample 

In order to investigate the relationship between MDO 
mass and nuclear radio luminosity, we have expanded 
our distance- limited (D < 19 Mpc) LLAGN sample as fol- 
lows. We include 73 additional Seyferts and LINERs (at 
D > 19 Mpc) from the hst of LLAGNs in H97a; these 
have also been observed at 2 cm and 150 mas resolution 
(Nagar et al. 2000, Nagar et al., in preparation). Also, 
we include an additional 31 galaxies which have reliable 
MDO mass estimates published in the literature. These 



include all galaxies in Richstone et al. (199S), Gebhardt 
3t al. (2000D, [Ferrarese fc Merritt (2000] ), and [Merritt fc 
Ferrarese (2001), and individual galaxies from Shields et 
al. (2000D and [Bower et al. (2000] ). Of the 200 galaxies, 154 
have all three of the following: an estimated MDO mass 
(either directly from dynamic measurements or estimated 
from (Tc), a core radio flux, and a bulge luminosity. The 
distribution of galaxy distance for these 154 galaxies is 
shown in Fig. |l|b. 

3. New Observations and Data Reduction 

3.1. VLA Observations and Data Reduction 

Seventy two LLAGNs at D < 19 Mpc were observed at 
2 cm (15 GHz) with the VLA. The bulk of the new obser- 
vations was carried out on 1998 March 3 and 1998 May 
21 while the VLA was in A-configuration ( Thompson et 
|al., 1980 ) and on 1998 June 21 while the array configura- 
tion was being changed from A- to B-configuration. Data 
taken during the latter period have more or less the same 
resolution as A-configuration data as the outermost an- 
tennas of the array had not yet been moved. Six LLAGNs 
(IC 239, IC 1727, NGC 428, NGC 660, NGC 1055, and 
NGC 1058) were observed on 1998 July 31 while the ar- 
ray was in B-configuration and two LLAGNs (NGC 3941 
and NGC 7177) were observed on 1998 September 5 and 
September 10, respectively, while the array was in A- 
configuration. Most of the data (all data taken before July 
1998) were calibrated and mapped by HF and preliminary 
results have been presented in Falcke et al. (1999). The re- 



maining data were calibrated and mapped by NN. 



In view of the non-simultaneity of the radio data used to 
derive the spectral index, we use a — —0.3 instead of —0.5 as 
the cutoff between a steep and flat spectrum. 
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Data were calibrated and mapped using AIPS, follow- 
ing the standard reduction procedures as outlined in the 
AIPS cookbook^. Observations of 3C 286 and 0404-^768 
(for which we used a 2 cm flux-density of 1.46 Jy) were 
used to set the flux-density scale at 2 cm. For sources 
stronger than about 3 mJy, we were able to iteratively self- 
calibrate the data so as to increase the signal-to-noise ra- 
tio in the final map. The typical root mean square (r.m.s.) 
noise in the final maps was 0.2 mJy and we take the detec- 
tion limit of the survey to be 1 mJy (i.e. 5 a). The resolu- 
tion of the final maps was ~ 0'.'I5 for the A-configuration 
maps and ^ 0'.'4 for the B-configuration maps. 

The 1(7 error in flux bootstrapping (i.e. setting the 
flux density scale relative to the flux calibrators) is ex- 
pected to be roughly 2.5%. Elevation dependent effects 
(atmospheric opacity and varying antenna performance 
with elevation) lead to additional errors when measuring 
the fluxes of some sources (see Sec. |4^ ). Such elevation de- 
pendent effects are expected to be small for most sources 
as the flux calibrator sources were observed at elevations 
45° to 50° and most sources were observed at elevations 
between 40° and 60°. We therefore did not correct for ele- 
vation dependent effects, but take them into account when 
estimating the error in the flux measurement. 

3.2. VLBA Observations and Data Reduction 

Nine of the sixteen nuclei listed in Table |^ have been pre- 
viously observed with the VLBA or VLBI (see column 5), 
and new VLBA observations of the remaining seven are 
presented here. We also observed two additional LLAGNs 
(with D > 19 Mpc) from the Palomar sample. The first, 
NGC 2655, has S^^^f > 2.7 mJy, but does not have a 
flat-spectrum radio core (Nagar et al., 2000). The other. 



NGC 3147 has both S^^,'^ > 2.7 mJy an d a flat-spectnmr 
radio core (Nagar et al. 
2001bD. 



2000, see also Ulvestad & Ho, 



These nine LLAGNs were observed at 6 cm (4.9 GHz) 
in a 10 hour VLBA run on April 1, 1999. Each source 
was observed at two different hour angles in order to ob- 
tain good (u, w)-coverage. At each hour angle, a cycle of 
6.5 min on source -I- 1 min on a nearby phase calibra- 
tor was repeated four times. The total integration time 
on each source was therefore 52 minutes. The "fringe 
finder" sources 4C 39.25 and 3C 345 were also briefly ob- 
served, and a 7 minute observation of J1642-I-6856 was 
used to check (and correct) the amplitude calibration of 
the VLBA. The weather was fair (and dry) at all VLBA 
sites, except for occasional wind gusts in excess of 35 mph 
at Kitt Peak. There were a few intermittent tape prob- 
lems at several antennas; all data points with tape weights 
less than 0.7 (on a scale of to 1) were flagged as bad. 
Data were calibrated using AIPS, following the standard 
procedures outlined in chapter 9 of the AIPS cookbook. 
Bad (li, v) data were deleted before the phase solutions of 



the phase-calibrator observations were transfered to the 
galaxy data. 

The sources were initially imaged using AIPS task 
IMAGR. For sources stronger than about 4 mJy, we were 
able to iteratively self-calibrate and image the data so 
as to increase the signal-to-noise in the final map. The 
self-calibration and imaging cycles were done using the 
DIFMAP package (^heperd, 1997|) , and the final self- 
calibrated (m, v) flies were then imaged using the AIPS 
task IMAGR. The peak flux-density of the source typically 
increased by a factor of 1.3 during the self-calibration pro- 
cess. Therefore, for sources weaker than 4 mJy, on which 
accurate self-calibration was not possible, we have mul- 
tiplied the peak detected flux-density by 1.3 as a crude 
attempt to correct for atmospheric decorrelation losses. 
The r.m.s. noise in the final uniformly weighted images is 
typically 0.15 mJy to 0.2 mJy, and the resolution between 
2 mas and 5 mas. 



4. Observational Results 

4.1. Results of the New VLA Observations 

The results of the new VLA observations are listed in 
Table || with columns explained in the footnotes. All 
six LLAGNs observed in B-conflguration (0'.'4 resolution) 
were undetected. All remaining LLAGNs have been ob- 
served in A-configuration or at a r esol ution equivalent to 
that of the A-conflguration (Sect. 3.1). Thus, the fluxes 
for all radio detected sources are measured at ~ 0'.'15 res- 
olution. The radio positions for the detected nuclei are 
limited by the positional accuracy of the phase calibra- 
tors, which is typically 2-10 mas, and by the accuracy 
of the Gaussian fit to the source brightness distribution, 
which depends on the signal to noise ratio of the source 
detection. The overall accuracy should typically be better 
than ~ 50 mas. We have compared the radio positions de- 
rived here with optical positions from Cotton, Condon, & 
Arbizzani (1999), which were measured from the digital 
sky survey with typical la accuracy l'.'5-2'.'5 in each of 
right ascension and declination. The results (column 6 of 
Table |^) show a good (< 2a) agreement in most cases. 

Only two nuclei in Table |, NGC 3628 and NGC 4486 
(M 87), show reliable extended structure in our 2 cm maps. 
The extended emission in the starburst galaxy NGC 3628 
is from a known star-forming region (Carral, Turner, & 



Ho, 1990| ) whi le that in NGC 4486 is from the well-known 
radio jet (e.g. Junor fc Birctta, 1995). The absence of ex- 



available online at www.nrao.edu 



tended emission in the other nuclei is not surprising as 
the high resolution may resolve it out; in addition, such 
emission is expected to be weak at the high frequency ob- 
served. For a few sources, a Gaussian fit (with peak fiux- 
density, major and minor axes as free parameters) to the 
nuclear radio source brightness distribution resulted in a 
Gaussian size slightly larger than the beam size so that the 
peak flux-density was slightly smaller than the total flux; 
however, the deconvolved source sizes were much smaller 
than half the beam size, so we consider these sources as 
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unresolved. Thus, most of the detected 2 cm nuclear radio 



sources are compact at the O'.'IS resolution (typically 15- 
25 pc) of our survey. The implied brightness temperatures 
for the 2 cm compact nuclei are typically T^ > lO^-^"^ ^ K. 

Non- simultaneous spectral indices for the 2 cm de- 
tected nuclei have been estimated by a comparison with 
radio data at other wavelengths from the literature (as 
indicated in column 12 of Table |^). While the actual 
value of the spectral index is uncertain given the reso- 
lution mismatch and the non-simultaneity of the observa- 
tions, we can determine whether the core is flat spectrum 
(a > —0.3; Si, cx i^"), or steep spectrum (a < —0.3), 
as noted in column 11 of the table. We are confident of 
the reality of the flat spectra obtained for the following 
two reasons. Since compact flat-spectrum radio sources 
are often variable, the use of non-simultaneous data at two 
frequencies can cause some intrinsically flat-spectrum ra- 
dio sources to be misclassified as steep-spectrum sources. 
However, since extended steep-spectrum radio sources are 
not variable, the use of non-simultaneous data at two fre- 
quencies should rarely cause intrinsically steep-spectrum 
radio sources to be misclassifled as flat-spectrum sources. 
Also, since the resolution is better at the higher frequency, 
resolution effects will tend to steepen the measured spec- 
trum if extended emission is present. Six of the 2 cm de- 
tected sources are not marked as flat-spectrum sources 
in the table. Of these, three show extended radio emis- 
sion from jets or lobes - NGC 4388 (|Falcke et al., 1998| ), 



LLAGNs was also observed in October 1996 (Nagar et al 
2000). The 2 cm emission from several other LLAGNs was 



NGC 4636 (^tanger fc Warwick, 1986^ , and NGC 4472 
( Ekers fc Kotanyi, 197^ ), and one shows extended emis- 
sion from star-formation - NGC 3628 (Carral, Turner, & 
Ho, 1910). This extended emission may dominate over any 
flat-spectrum nuclear component within the beam of our 
radio maps. For the two remaining sources, NGC 4293 
and NGC 4419, VLA archive data at 6 cm (0'.'4 resolution; 
peak flux-densities 4.1 mJy/beam and 11.8 mJy/beam, re- 
spectively) and 20 cm (1'.'2 resolution; peak flux-densities 
13.0 mJy/beam and 28.3 mJy/beam, respectively) show 
that both have steep-spectrum radio cores which are more 
or less unresolved at both frequencies. 



4.2. Variability of the VLA-Observed Radio Emission 

Except for a few well-studied cases, e.g. NGC 3031 [M 81] 



( Ho et al., 1999| ; |Bietenholz et al., 20001 ) and NGC 5548 



( Wrobel, 2000 ), little is known about radio variability in 
Seyferts and LLAGNs. In the course of our observations, 
we have observed several LLAGNs at multiple epochs 
using the same telescope (VLA), telescope conflguration 
(A conflguration), wavelength (2 cm), and data reduc- 
tion methods ( Nagar et al., 2000| , 2001, Nagar et al., in 



preparation, this work). This offers us a unique opportu- 
nity to study radio variability in a sample of LLAGNs. 
Specifically, the 2 cm emission from all sixteen LLAGNs 
listed in Table |l| was observed at two epochs - March, May 



also observed at these epochs. In addition, the 3.6 cm emis- 
sion from the sixteen LLAGNs in Table ^ was observed in 
September 1999 ( [Nagar et al., 200lD and December 2000 
(Nagar et al., in preparation). 

The errors in the measured fluxes come from five signif- 
icant factors (see the VLA Observational Status Summary 
for details^ (a) the accuracy of flux bootstrapping from 
a VLA 'flux' calibrator. The 3.6 cm fluxes at both epochs 
were bootstrapped from observations of 3C 286. The 
2 cm fluxes were bootstrapped from observations of either 
3C 286 and 0404-^768 (for which we used a 2 cm flux- 
density of 1.46 Jy at all epochs). All VLA-recommended 
flux bootstrapping procedures were followed and the (la) 
errors due to this factor are estimated to be 1% at 3.6 cm 
and 2.5% at 2 cm. (b) variation of antenna gain with ele- 
vation. The gain of the VLA antennas at 2 cm and 3.6 cm 
is maximum at elevation Eg = 45°-50° and decreases (i.e. 
decreasing the estimate of the observed flux) as the ele- 
vation increases or decreases from Eq. At both epochs, all 
3.6 cm observations were made at elevations 35°-70° with 
the flux calibrators observed at elevation ^40°. The an- 
tenna gain variation over this elevation range is expected 
to be <0.3% and we therefore ignored this effect at 3.6 cm. 
The antenna gain variation at 2 cm is larger; it is less than 
1.5% over the elevation range 35° to 65°, but increases to 
^--^4% at elevation 80°. The 2 cm flux calibrators were ob- 
served at elevations between 40° and 55° in all cases. The 
variation of gain with elevation was corrected for only in 
the reduction of the September 1999 data. For the other 
2 cm data we took this effect into account when estimat- 
ing the error of the flux determination, (c) atmospheric 
opacity. The atmospheric opacity at the VLA is typically 
0.02 at 2 cm and <0.01 at 3.6 cm. At these opacities, 
atmospheric absorption will result in a source having an 
observed flux which is < 0.5% larger when observed at el- 
evation 70° than at elevation 45°. We did not measure 
the atmospheric opacity during any of the observing runs 
except that of December 2000. From the weather condi- 
tions during the observing runs it is reasonable to assume 
the opacities were not significantly higher than the typical 
values quoted above. We did not correct for atmospheric 
opacity effects in the 2 cm and 3.6 cm data. Therefore, 
for sources observed at elevations below 45° we use the 
typical opacities above to estimate the additional error in 
the measured flux; It is notable that in the elevation range 
45°-70° this effect works opposite to the effect discussed 
in (b), thus decreasing the overall error, (d) the r.m.s. 
noise in the final maps, which was typically 0.1 mJy and 
0.2 mJy at 3.6 cm and 2 cm, respectively, (e) uncertainties 
intrinsic to the self-calibration process itself. 
In summary, factor (a) usually dominates so that the la 
fiux calibration errors are typically about 1% and 2.5% for 
3.6 cm and 2 cm observations at the VLA, respectively. 
While we have taken some effort to estimate the flux er- 



or June 1998 (this work), and September 1999 (Nagar et 

al., 2001). The 2 cm emission from some of these sixteen www.nrao.edu 
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rors, other factors, e.g. variability in the flux calibrators, 
may cause the errors to be higher than estimated here, 
especially at 2 cm. 

The inter-year variability for the LLAGNs with multi- 
epoch observations is shown in Fig. |^. We have plotted 
<S> versus AS/<S>, where AS is the r.m.s. flux varia- 
tion and <S> the average flux (for nuclei with only two 
epochs of measurements, AS represents the difference be- 
tween the fluxes). The grey shading delineates the region 
within which observed variability can be considered reli- 
able. The lower cutoff to the shaded region is at 5 mJy. 
Sources with flux lower than this could not be accurately 
self-calibrated, which results in larger than usual errors in 
their flux determination. The left cutoff results from un- 
certainty in the flux-calibration. This uncertainty coupled 
with the r.m.s. noise in the radio maps could lead to values 
of AS/<S> which are ~ 10%, even if the source did not 
vary. As an alternate illustration of the flux variabilitv we 



tor positions have now been determined to an accuracy of 
1 mas or better (Bcasley et al., 2002) and the new posi- 
tions are available on the web|^ We therefore updated the 
positions of the sources to reflect the new accurate phase 
calibrator positions. The other factor contributing to the 
position uncertainty is the accuracy of the Gaussian fit to 
the source, which should typically be better than 2 mas. 
Thus the overall accuracy of the positions listed in Table || 
should be about 2 mas or better. The implied brightness 
temperatures were calculated using the formula given in 
Falcke et al. (2000| ); the results are in the range 10^'* to 
10^-^ K. Since most of the sources are unresolved, these 
values are lower limits to the true brightness temperatures. 

Among the newly observed sources, the two with the 
highest peak flux-density at 6 cm, NGC 4374 (M 84) 
and NGC 4552 (M 89), clearly show extended structures 
(Fig. ^) suggestive of parsec-scale "jets." T he pc-scale ex- 



tension in NGC 4374 (detected earlier by Wrobel et al 



have p lotted the 2 cm 'Ught curves' for five of the nuclei 1996D is aligned with the lar ger kpc-scale FR-I radio lobes 



in Fig. ||. In Fig. || the data for the first three epochs are 
from our 2 cm observations. The data for the last epoch 
shown (December 2000) has been obtained by scaling the 
previous epoch's 2 cm flux by the ratio of the 3.6 cm fluxes 
which were determined at both epochs. 

3an bo soon in the Figures, we find significant vari 



As 



ability l in almost half the sources . While the 2 cm flux val - 
ues potentially suffer from significant flux calibration and 
elevation-dependent errors, the 3.6 cm flux-calibration is 
quite reliable. The result is made more cr edible by earlier 



(e.g. [Laing fc Bridle, 1987| ), while the twin pc-scale ex- 
tensions in NGC 4552 are, to our knowledge, the first 
time such structure has been detected in this galaxy. As 
pointed out in Nagar et al. (2001), of the sixteen objects in 
Table |l| the five with the highest 6 cm peak flux-densities 
at mas resolution - NGC 3031 (M 81; Bietenholz et al. 
"2000| ), NGC 4278 ( [Jones et al., 1984|; ^^e et al., 2000| ), 
NGC 4 486 (M 87; [Junor fc Biretta, 1995| ), NGC 4374 
(M 84; [Wrobel et al., 1996], this work), and NGC 4552 



observations of variabi lity in NGC 3031 ( Ho et al., 199S 



Bietenholz et al, 200C| ) and NGC 4486 ( [Junor fc Biretta, 



1995|). [The large periods between the different epochs do 
not allow for detailed investigation of the variation, and at 
this point we are unable to find any clear correlation be- 
tween the amount of variability and the other properties 
of the LLAGNs. Future higher time-resolution, multifre- 
quency data will be useful in this context. 

4.3. Results of the New VLB A Observations 

All 7 objects from Table |l] which were newly observed 
with the VLBA at 6 cm were clearly detected in initial 
maps (i.e. without any form of self-calibration). Of the two 
other (D > 19 Mpc) LLAGNs observed, NGC 3147 was 
detected, but NGC 2655 was not detected at a lOcr limit of 
0.76 mJy. The non-detection of NGC 2655 is not surpris- 
ing as this object has a steep-spectrum radio core at arc- 
sec resolution (Nagar et al., 200C). Of the detected nuclei, 
mas-scale radio cores were known to exist in NGC 4374 
and NGC 4552 (e.g. [Jones, Terzian, fc Sramek, 198l|). The 



(M89; this work) - all have one- or two-sided pc-scale ex- 
tensions, suggestive of jets. For these sources the peak flux- 
density of the extended emission is typically 10%-25% of 
the peak flux-density of the core. The signal-to-noise ratios 
of the core detections for the other eleven sources listed in 
Table range between 10 and 50. The relative weakness 
of most of these sources does not allow very good self- 
calibration. This weakness, coupled with the low signal- 
to-noise expected for any extended component (10%-25% 
of a 10 fT-50 (T detection) may provide an explanation for 
the non-detection of weak extended structure in most of 
these sources. In fact, very de ep VLBA observations of 
NGC 4258 reveal a sub-pc jet ([Hcrrnstcin et al., 1997|). 



results of the new VLBA observations are listed in Table |, 
with columns explained in the footnotes. The source po- 
sitions are referenced to the positions of their respective 
phase-calibrators. All phase calibrators were drawn from 
the VLBA calibrator survey and the positions we used 
at the time of observation were expected to be accurate 
to ±10 mas in each of R.A. and Dec. The phase calibra- 



5. Detection Rates of Radio Cores in the 
Distance-Limited Sample 

5.1. Parsec-Scale Radio Cores (VLBA Detections) 

The new VLBA results conflrm that all LLAGNs at 
D < 19 Mpc with a 2 cm compact flat-spectrum core (as 
determined from VLA observations), i.e. 17% ±4% of 
all LLAGNs at D < 19 Mpc, have mas-scale radio cores 
with brightness-temperature ^ 10^ K. It is notable that 
this sample of LLAGNs with ultracompact radio cores is 
almost equally divided between Seyferts (6 of 16) and 
LINERs (8 of 16, see Table |l|). Nuclear starbursts can 
have a maximum brightness-temperature of ~ 10^~^ K 
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( |Condon et al., 1991 ) while the most luminous known ra- 
dio supernova remnants (e.g. Colina et al., 2001) would 
have brightness temperatures < 10*^ K even if they were 
< 1 pc in extent. As argued in Falcke et al. (2000), if 
the core radio emission is attributed to thermal processes, 
the predicted soft X-ray luminosities of LLAGNs would be 
at least two orders of magnitude higher than observed by 
ASCA ([Terashima, Ho, fc Ptak, 2000|) and Chandra po et| 



al., 20C|i|). Fi nally, as discussed in[Ulvestad fc Ho (2001a| ), 
single SNRs (polina et al., 2001] ) or a collection of SNRs 
( Neff fc Ulvestad, 200C| ) would have radio spectral indices 
a ~ —0.7 to —0.4 rather than the values a ~ —0.2 to 



0.2 seen in the VLBA-detected LLAGNs ( Nagar et al ^ 
200l|)~frhus, the only currently accepted paradigm which 
may account for the mas radio cores is accretion onto a 
supermassive black hole. In this case, the mas-scale radio 
emission is likely to be either e mission from the accre- 
tion inflow (Narayan et al., 1998) or synchrotron emission 
from the base of the radio jet launched by the accreting su- 



permassive black hole ( Falcke fc Biermann, 1999 ; Zensus, 
19970 

Fifteen of the sixteen definite AGNs (Table ^ are ei- 
ther (a) type 2 LLAGNs in massive (Mb (total) < -20) 
ellipticals (four galaxies), or (b) type 1 LLAGNs (two el- 
lipticals and nine non-ellipticals). The remaining galaxy, 
NGC 5866, is a transition nucleus in an SO galaxy with an 
ambiguous detection of broad Ha (H97c); the ambiguity 
is not unexpected as the dilution by H II regions in tran- 
sition nuclei makes it more difficult to determine the pres- 
ence of broad Ha than in "pure" LINERs and Seyferts (see 
H97c). NGC 5866 may, therefore, belong with the other 
nine type 1 non-ellipticals. 



5.2. 100-Parsec-Scale Radio Cores (VLA Detections) 

We have completed observations of 95 of the 96 LLAGNs 
in the distance-limited sample with the VLA; only IC 356, 
a transition nucleus, remains unobserved by us. IC 356 can 
be considered as a non-detection for our purposes since a 
significantly lower resolution (1') fi ux measurement gi ves 
an upper limit of 2 mJy at 2.8 cm ( Niklas et al., 1995 ). 

When all 96 LLAGNs at D < 19 Mpc are considered to- 
gether, our VLA observations have detected 9 of 23 (39%) 
low-luminosity Seyferts^, 16 of 37 (43%) LINERs, and 6 
of 35 (17%) transition nuclei, at a resolution of ~0'.'15 
and above a fiux limit of ^ 1 mJy (Nagar et al. 2000; 
this work) . Alternately, we can state that radio cores with 
luminosity P^^cm > 10^° W Hz^^ are found in 5 of 23 low- 
luminosity Seyferts, 8 of 37 LINERs, and 3 of 35 transition 
nuclei in the distance-limited sample. The radio luminosi- 
ties of the detected 2 cm cores lie between 10^^ and 10^^ 
W Hz~^ (e.g. Fig. H), similar to the luminosities seen in 
'normal' E/SO galaxies (Sadler et al. 1989). It is notable, 



however, that a significant fraction of the detected 2 cm 
compact cores are in spiral galaxies (Fig. ||a). 



6. Correlations between Radio Power and Optical 
Emission-Line Properties in the 
Distance-Limited Sample 

Emission-line fiuxes and linewidths for the galaxies in the 
Palomar sample were typically measured over a 4" x 2" 
aperture centered on the galaxy nucleus (H97a). In this 
section we study the relationships between these nuclear 
emission-line properties and the core radio power mea- 
sured at 150 mas resolution for LLAGNs in the distance- 
limited (D < 19 Mpc) sample. 

6.1. Emission-Line Luminosity 

Correlations between the radio power and emission-line lu- 
minosity of active galaxies are well known, and are found 
in both Seyfert (do Bruyn fc Wilson, 1978 ) and radio 
galaxies ( |Baum fc Hcckman, 1989| ). The [O III] A5007 line 
is usually used in assessing these correlations, but this 
line is weak in LINERs and we use instead [O I] A6300, 
which also has the advantage of being uncontaminated by 
emission from HII regions. In our D < 19 Mpc LLAGN 
sample the 2 cm radio power appears approximately pro- 
portional to the nuclear [O I] A6300 luminosity (Fig. ||) , 
though there is a range of about two orders in magni- 
tude of emission-line luminosity for a given radio power. 
Testing the statistical significance of the correlation is not 
straightforward as there are some uncertain values (non- 
photometric data) of the [O I] luminosities. For simplic- 
ity, we treat the non-photometric flux densities as photo- 
metric measurements when using the bivariate tests from 



NGC 185 is considered a non-detection since the radio 



source detec ted in this galaxy is probably extra-nuclear ( Nagai 
et al, 20001). 



the A SURV statistical software package (Lavalley et al. 
1992 , hereafter Asurv). This approximation should not 
bias the results significantly as the uncertainties of the 
non-photometric emission-line fiuxes are expected to be 
a factor of only ±2 (H97a), which is small compared to 
the 3 orders of magnitude range of the [O I] luminosity. 
When this is done, the Asurv tests indicate that the [O I] 
luminosity is correlated with the 2 cm radio power at 
the 98.5%-99.9% significance level for the thirty-one 2 cm 
detected LLAGNs at D < 19 Mpc, at the 99.9% signifi- 
cance level for all Seyferts and LINERs at D < 19 Mpc, 
and at the 99.99% significance level for aU 95 LLAGNs 
at D < 19 Mpc which have been observed at 2 cm (see 
Table ^ . The correlation among the radio-detected nuclei 
is not just a result of the elliptical galaxy nuclei having 
higher radio powers and emission-line luminosities than 
the non-ellipticals; among the 24 non-elliptical radio de- 
tections, the radio power and [O I] luminosity are corre- 
lated at the 99.5% significance (Table ||). 

Interestingly, if we use the Ha+[N II] AA6548,6583 
luminosity in place of the [O I] luminosity the data ap- 
pear more scattered, and the significance of the correla- 
tion drops to 97%-98% for all radio detected LLAGNs and 
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89%-93% for the radio-detected non-elliptical LLAGNs 
(Table |). 

6.2. Emission-Line Width 

The only emission line for which the full width at half max- 
imum (FWHM) is reported in the Palomar survey results 
is [N II] A6583 (H97a). We therefore use this line for assess- 
ing the correlations between linewidths in LLAGNs and 
other parameters. For all radio detected nuclei considered 
together, aU three of the [Nil] A6583 FWHM linewidth, 
[N II] A6583 luminosity, and 2 cm radio power are corre- 
lated, with the significance of the correlations as follows: 
linewidth and radio power 99.9%, linewidth and [Nil] lu- 
minosity 99-99.5%, and [Nil] luminosity and radio power 
(97-99%; T able 3). The sample of low- luminosity Seyferts 
studied by Ulvestad fc Ho (2001a ) also follows the above 



detected nuclei with [Nil] linewidths in excess of AV^r^j 
(Figs. I and 0) are NGC 4278, NGC 4374 (both ellipticals 
with pc- to kpc-scale radio jets), and NGC 4772 which has 
a pc-scale radio core (Sect. 



three correlations at a high significance level. The [Nil] 
linewidths of LINERs and transition nuclei are generally 
larger for the earlier type galaxies (the Seyferts do not 
show such a clear variation with morphological type). 
Also, the same early type galaxies tend to have higher 
radio powers though not necessarily higher [Nil] lumi- 
nosities. Interestingly, the [Nil] linewidth does not show 
any clear dependence on the bulge or total magnitude of 
the galaxy among the LLAGNs with D < 19 Mpc. 

That the galaxy rotation contributes to the [Nil] 
linewidth is clear: the [Nil] linewidth is correlated with 
the corrected (for galaxy inclination) rotational velocity 
of the host galaxy (AV^J^t) for all LINERs (99.7% signifi- 
cance) and transition nuclei (97%-98% significance) with 
D < 19 Mpc (Fig. I and Table |). The significance of the 
above correlation for only the radio detected LINERs is 
smaller (82%-95%) perhaps due to the small number of ob- 
jects. More interestingly, there is evidence that the nuclear 
radio power is related to the ratio FWHM([N II])/AV;:ot, 
which increases with increasing radio power (Fig. [7|). This 
correlation has a significance of 99.9% if all radio de- 
tections are considered, and 98%-99% if only the non- 
elliptical radio detections are considered (Table ||). That 
this correlation may involve a direct physical significance 
is supported by the lack of a clear correlation between 
[Nil] width and bulge luminosity among the LLAGNs 
with D < 19 Mpc. One immediate explanation for the cor- 
relation is that the radio source is responsible for "stir- 
ring" the [Nil] e mitting gas as is probably the case wit h 
powerful Seyferts ( |Wilson fc Willis, 1980| ; [Whittle, 1992b| ). 
The converse, however, is not true. That is, a high ratio of 
FWHM([N II])/AVJot does not necessarily imply a strong 
radio source: notice how the radio non-detections (Fig. |^) 
show a large variation in the y-axis. This is one of the few 
suggestions that the radio non-detected nuclei are not sim- 
ply lower radio-power versions of the radio detected nuclei. 
In both Figs. |and 0, the radio-detected nuclei NGC 4143, 
NGC 4168, NGC 4486, NGC 4552, and NGC 5866 are not 
plotted as these nuclei do not have values for the galaxy 
rotational velocity listed in H97a. The only three radio 



7. Comparison of the Radio Detected and 
Non-Detected LLAGNs in the 
Distance-Limited Sample 

The sixteen LLAGNs observed at mas resolution and dis- 



cussed in Sect. 5.1 represent a flux-limited sample drawn 
from the 2 cm VLA-detected LLAGNs. It is likely that 
some fraction of the VLA-detected radio cores which have 
not yet been observed with the VLBA, and also some frac- 
tion of the VLA-non-detected LLAGNs, also have high- 
brightness temperature cores (and thus accreting black 
holes). These issues, crucial for determining the true in- 
cidence and properties of accreting black holes in nearby 
galaxies, are explored in the following sections. The fluxes 
and detection rates quoted in the following sections are 
those from the VLA survey. 

7.1. LINERs and Low-luminosity Seyferts 

The type 2 nuclei (broad Ha not detected) have a rather 
low VLA radio core detection rate: only 1 of 11 Seyfert 
2 nuclei and only 8 of 27 LINER 2 nuclei are detected. 
The type 1 nuclei (broad Ha present) have a significantly 
higher detection rate: 8 of 10 LINER 1.9s, and 8 of 12 
Seyfert Is (1.0s to 1.9s) are detected at 2 cm. The differ- 
ence between the detection rates of type Is and type 2s 
increases when one considers the type classification more 
carefully. Of the 22 type 1 nuclei in the distance-limited 



sample, a definite detection (sec Ho et al., 1997c) of broad 
Ha has been made in 17; the remaining five have probable 
detections of broad Ha. Also, 6 of the 38 LINER type 2 
and Seyfert type 2 nuclei in the distance-limited sample 
have ambiguous ( Ho et al., 1997c ) detections of broad Ha. 
Of the definite type Is, 13 of 17 (76%) are detected in 
the radio. The only four non-detected definite type Is 
are NGC 2681, NGC 4051, NGC 4395, and NGC 4639. 
Excluding type 2s with an ambiguous detection of broad 
Ha, 7 of 32 (22%) type 2 LINERs and Seyferts are de- 
tected in the radio. 

It is known that the radio powers of AGN are cor- 
related with various parameters of the host galaxy, such 
as the bulge luminosity (e.g. Nelson fc Whittle, 1996 ) and 
possibly the nuclear luminosity at optical/UV wavelength s 
(e.g. |de Bruyn fc Wilson, 1978|; |Baum fc Heckman, 1989|). 



It is, therefore, important to see if the preferentially higher 
detection rate of the type 1 nuclei over the type 2s is influ- 
enced by these effects. The single type 2 Seyfert detected 
in the radio, NGC 4472, is a massive (Mb (total) < -20) el- 
liptical. Of the eight type 2 LINERs detected in the radio, 
two (NGC 4374 and NGC 4486) are massive eUipticals, 
one (NGC 2841) has an ambiguous identification of broad 
Ha(H97c), four have very weak (~ 3-7 cr) radio detections 
and the last (NGC 4736, 1.9 mJy or ~ 10 cr, distance 
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4 Mpc) is unusually nearby and would not be detected 
at the median distance of the sample. Turning to the 
type Is, the two non-detected LINER 1.9s are NGC 2681 
and NGC 4438. The former has among the lowest bulge 
luminosity, emission-line luminosities and [Nil] FWHM 



Further support for this idea comes from an apparent pro- 
portionality between core radio power and luminosity of 
the broad Ha component (Fig. ||a) for all (20) type 1.9s 
in the Palomar sample for which broad Ha luminosities 



were measured under photometric conditions (Ho ct al 



in the LINER 1.9 sample, and the broad Ha detection 1997c). Twelve of the plotted nuclei are in the distance- 



is doubtful in the latter (H97c). The four non-detected 
Seyfert type Is are NGC 3982, NGC 4051, NGC 4395, and 
NGC 4639. The broad Ha detection in NGC 3982 is doubt- 
ful (H97c) and the narrow Ha luminosity of NGC 4051 
places it exactly on the operational cutoff between low- 
luminosity Seyferts and 'classical' Seyferts, so it may have 
a different type of central engine from LLAGNs. A sub- 
mJy mas-scale nuclear radio source with brightness tem- 
perature > 2 X 10^, potentially an accreting massive black 



hole, has been found in NGC 4395 ( Wrobel, Fassnacht, fc 
Ho, 20^. Such a weak radio core is expected in view of 
the extremely low bulge luminosity (MB(total) = —10.4) of 
this galaxy. In summary, the radio-detected type 2s tend 
to be strong radio sources in high bulge luminosity galax- 
ies or very weak radio sources, and of the six type 1 nuclei 
not detected in the radio, two have a doubtful detection 
of broad Ha and two have unusually low bulge luminosi- 
ties which may be the reason for their radio non-detection. 
The preference for radio detected LLAGNs to be in type 1 
nuclei is the only highly significant difference between the 
2 cm detected and non-detected LINERs and Seyferts. 
Thus, it is most likely that the difference in radio power 
is indeed related to the AGN type. 

In this context it is relevant to look at the detection 
rate of the elliptical LLAGNs at D < 19 Mpc. AU four mas- 
sive (Mb (total) < —20) type 2 ellipticals (one of which 
has a transition nucleus) and all three type 1 ellipticals 
(with MB(total) values of -19.1, -18.9, and -20.7) have 
detected radio cores. Only three ellipticals (all type 2) re- 
main undetected in the radio: NGC 185, NGC 3379, and 
NGC 4494. NGC 185 is a nearby very low bulge mass 
(Mb (total) = —14.9) elliptical and is thus expected to have 
a very low luminosity radio core; NGC 33 79 (Mb (total) — 
— 19.4) may actually be an SO galaxy ( Statler, 2001 ); 
NGC 4494 (MB(total) = -19.4) has a very low emission- 
hne luminosity (log [Lhq/W] < 30.54; H97a). As noted 
above, the radio detection of the massive ellipticals would 
be expected if the radio power scales with the bulge lumi- 
nosity in LLAGNs (Sect. B.2). Among less massive ellip- 
ticals it may be that the type 1 nuclei are preferentially 
detected similar to the case for non-ellipticals. 

Thus the demography of LLAGNs with VLA-detected 
radio cores is similar to that of the subset with VLBA- 
detected radio cores, i.e. radio cores are preferentially seen 
in either massive type 2 ellipticals or in type 1 LLAGNs. It 
is likely that we have detected 2 cm radio cores and broad 
Ha lines in only the more luminous LLAGNs, since the 
[O I] A6300 luminosities of the Seyfert and LINER type 1 
nuclei are higher than those of the type 2 nuclei at the 



limited sample; radio data for the remaining nuclei were 
taken from Nagar et al. (in preparation). The proportion- 
ality is unlikely to be due to distance as it is also seen in 
the equivalent plot of fluxes (Fig. ||b) . 

Other marginally significant differences between the 
radio-detected and radio-non-detected nuclei include the 
following: as compared to the non-detected LINERs and 
Seyferts, the 2 cm detected LINERs and Seyferts have 
higher [O I] luminosities (~ 93% significance level; see 
Fig. H), earlier morphological types (Fig. ||a), higher bulge 
luminosities (Fig. |lo| ), and reside in regions with a higher 
galaxy density (90% -95% signifi cance level; Fig. ||b). The 



galaxy density, from [TuUy (19881) and listed in H97a, is the 
density of all galaxies brighter than Mb = — 16 mag in the 
vicinity of the object of interest. In other words, 50% of all 
2 cm non-detected Seyferts and LINERs reside in galax- 
ies which are located in regions of local galaxy density 
<0.5Mpc~'^, as compared to only 16% for 2 cm detected 
Seyferts and LINERs. The difference in the galaxy density 
is primarily a result of galaxies with more massive bulges 
being in regions of higher galaxy density: the three most 
massive ellipticals in the sample (all radio-detected) are at 
the centers of clusters (pgai > 3) while the four galaxies 
with the least massive bulges are in regions with pgai < 1. 
When these very high and very low bulge mass galaxies 
are excluded, the difference in galaxy densities between 
the detected and non-detected galaxies is no longer signif- 
icant. A comparison of the ratio of radio to [O I] luminosity 
(Fig. ^ for radio detected and non-detected LLAGNs re- 
veals that the ellipticals in the sample tend to have higher 
Lrad/L[Oi] ratios than the non-ellipticals for both type 1 
and type 2 nuclei. The highest ratios are in the most mas- 
sive type 2 ellipticals. Our radio detection limits are not 
low enough to reveal any clear difference in this ratio be- 
tween the radio detected and non-detected nuclei. 

Most of the above factors are consistent with the idea 
that our radio observations have detected the more lu- 
minous LINERs and low-luminosity Seyferts, and that 
deeper radio observations of the sample will reveal a higher 
incidence of compact radio cores. 

7.2. Transition Nuclei 

The only|^ transition nucleus (out of 64 in the Palomar 
sample) which may have detected broad Ha emission is 
NGC 2985 (at D= 22.4 Mpc); this detection is doubtful 
(H97c) . Thus it is most appropriate to compare transition 
nuclei to other type 2 LLAGNs in the sample. Although 



99.9% significance level. Ulvestad & Ho (2001a) arrived 



at a similar conclusion in their analysis of the radio pow- 
ers of ~ 45 Seyferts drawn from the Palomar sample. 



^ NGC 1161, tentatively identified as a transition nucleus 
(H97a) and listed in H97c as having a definit e broad Ha line 
is not in the Palomar sample (Ho et al., 1995). 
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the del action rate of Seyferts and LINERs of both types 
is much higher than that of transition nuclei, the detec- 
tion rate for type 2 LINERs and type 2 Seyferts (10/38 
or 26%) is similar to that for all transition nuclei (6/35 or 
~ 17%). Asurv tests indicate that the 2 cm radio powers 
and the [O I] A6300 luminosities of the transition nuclei 
are not significantly different from those of the type 2 
LINERs and type 2 Seyferts. One way to test for a differ- 
ence between the central engines of transition nuclei and 
other type 2 LLAGNs is to look for a systematic differ- 
ence in the ratio of radio to [O I] A6300 luminosities, since 
[O I] A6300 is expected to trace the 'LINER' component 
rather than the star formation component. Our data do 
not show any clear difference (Fig. |o|) , but in view of the 
large number of upper limits to j^Prad/[OI] A6300 deeper 
radio imaging i s req uired to address this issue. Not sur- 
prisingly (Sect. 7.1), if we compare the transition nuclei 
to both type 1 and type 2 LINERs and Seyferts, then 
Asurv tests indicate that the 2 cm radio powers of the 
transition nuclei are lower at the 99.5% significance level. 
The [O I] A6300 luminosities of the transition nuclei are 
also lower than those of the Seyferts and LINERs in the 
sample at the 97%-99% significance level. 

Are there any clear differences between the 6 detected 
and the 29 non-detected transition nuclei? NGC 4552 is a 
massive elliptical and NGC 4419 and NGC 5866 have am- 
biguous identifications of broad lines (II97c) which may 
(see Sect. 7.1) be related to their radio core detections. 
Another likely factor is differences in the emission-line 
ratio s (Fig. The ratios plott ed in Fig. ^ are used 
(e.g. Veilleux & Osterbrock, 1987) to distinguish between 
Seyferts, LINERs and Hll-region type nuclei. The dot- 
ted lines show the spectral classification cutoffs used by 
H97a. Clearly the transition nuclei detected at 2 cm (filled 
squares) are preferentially found among transition nuclei 
which lie close to the LINERs and Seyferts in the plots 
(while the [NIIJ/Ha and [SIIJ/Ha ratios were not used 
for spectral classification it is clear that transition nuclei 
do not attain as high values of these ratios as LINERs and 
Seyferts). Two of the detected transition nuclei, NGC 3627 
and NGC 3628, are an interacting pair, with the activ- 
ity classification of NGC 3627 bordering on a Seyfert 2.0 
(II97a). The classification of the elliptical NGC 4552 as 
a transition nucleus is uncertain as its II/3 flux is uncer- 
tain at the 30%-50% level (II97a). The median morpho- 
logical type is T™°'^'''" = 3 for the undetected transition 
nuclei, while the detected transition nuclei have T val- 
ues of —5, —1, 1, 3, 3 and 3. Finally, there are no sig- 
nificant differences in the [O I] luminosities and the lo- 
cal galaxy density between the detected and non-detected 
transitions. Therefore, it appears that the strongest factor 
which determines the probability of detecting a radio core 
in a non-elliptical transition nucleus is the proximity of its 
emission-line ratios to "pure" LINERs and Seyferts. This 
gives support to the view that at least some transition 
nuclei have both "pure" LINER/Seyfert and HII region 
components, with the detectability of the radio core de- 



et al. (2000) have detected arcsec-scale, 3.6 cm cores with 



a flat radio spectrum, in 5 of 25 transition nuclei from 
the Palomar spectroscopic sample (two of their five cores, 
NGC 3627 and NGC 4552, are also found to be compact 
and flat spectrum by us). They find that the detected tran- 
sition nuclei are preferentially in early-type galaxies; our 
results are consistent with this. 

7.3. H II Nuclei 

While we have not investigated HII nuclei from the 
Palomar sample in the radio, it is instructive to check in 
what ways galaxies with HII region nuclei may differ from 
those with LLAGNs. As shown by |Ho et al. (1997b| ), the 



H II nuclei in the Palomar sample are in later-type hosts 
as compared to the LLAGNs. This is also true for the 
subset of the Palomar sample at D < 19 Mpc. Further, the 
median absolute magnitude in the B-band of the bulge for 
the HII region galaxies at D < 19 Mpc (Mg"^^''^" (bulge) = 
— 16.83) is significantly fainter than the equivalent me- 
dian value for the low-luminosity Seyferts and LINERs 
at D<19 Mpc (Mg<='^'^" (bulge) = -18.65). The median 
values have been obtained from the values of Mb (bulge) 
listed in H97a. The differences between the D < 19 Mpc 
galaxies in the Palomar sample with H II nuclei and those 
with LLAGNs are thus deeper than just their nuclear spec- 
tral classification. Thus, while the HII region nuclei may 
also harbor accreting black holes, these will form a dif- 
ferent population from those in LINERs, low-luminosity 
Seyferts, and transition nuclei. 

8. The Influence of the Galaxy Bulge on the 
Properties of the Distance- Limited Sample 

The bulge luminosity of the host galaxy is known to 
strongly influence the probability of finding a nuclear radio 
source in radio galaxies (Auriemma et al., 1977), and the 
power of the nuclear radio emission in early-type galax- 
ies (Sadler et al., 1989). Further, for radio and Seyfert 
galaxies, the nuclear centimeter radio luminosity at arcsec- 
scales is correlated with the bulge luminosity over six or- 
ders of radio power (Nelson & Whittle, 199(:). The bulge 
luminosity is also known to affect the nuclear emission- 



termined by the strength of the former component. Filho 



line luminosity in Seyfert galaxies (Whittle, 1992b), radio- 
galaxies (Zirbel & Baum, 1995), and even "normal" early- 
type galaxies (Sadler et al., 1989). In this section, we test 
whether the above relationships hold at the lower radio 
and emission-line luminosities of our sample of LLAGNs. 

8.1. Emission-Line Luminosity 

Among the radio detected nuclei, the [O I] A6300 lumi- 
nosity is not clearly related to the galaxy bulge magni- 
tude (Fig. |l|). Also, for -17 > MB(bulge) > -20 the 
radio detected nuclei and non-detected nuclei are fairly 
well mixed at any given bulge magnitude. It is only at 
the lowest and highest bulge magnitudes that all nuclei 
are undetected and detected in the radio, respectively. 
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Also shown in the figure are the low-luminosity extrap- 
olations of the relationship between Lpi] and Mb (bulge) 
for 'classical' Seyferts (dotted line; for these nuclei the 
[O III] A5007 luminosity from [Whittle (1992a| ) was con- 
verted to the [O I] luminosity using a factor of 0.04, typ- 
ical of Seyferts ( |Veilleux fc Osterbrock, 1987 )) and an 
indicative line for FR I radio galaxies (dashed line; the 
total Ha+[N II] AA6548,6583 luminosity from [Zirbel fc 



the plot of L[oi] vs Mb (bulge) given in Fig. |T^, Fig. [l^ 
shows that for a given galaxy total or bulge luminosity, 
the radio detected and radio non-detected galaxies have 
similar nuclear emission-line luminosities. It is only at the 
lowest and highest bulge luminosities that all galaxies are 
undetected and detected, respectively, in our radio survey. 



Baum [(1995| ) was converted to an [01] luminosity us- 8.2. Core Radio Power 
ing a factor of 0.08, typical of the LINERs in our sam- 
ple), ^irbcl fc Baum (1995 ) list total emission-line lu- 
minosities for radio galaxies; typically more than half of 
t his luminosity is from the nuclear (<2.5 kpc) region 
( Baum fc Hcckman, 1989 ). It is notable that, while several 
Seyfert galaxies fall close to the low-luminosity extrapo- 
lation of 'classical' Seyferts, many lie much further away, 
the farthest being NGC 4168, NGC 4472 (both ellipticals), 
NGC 4565, and NGC 4725. The elliptical LINER with a 
high [O I] luminosity for its bulge magnitude is NGC 4278. 
The radio-detected Seyferts are indistinguishable from the 
other radio-detected LLAGNs and, when considered to- 
gether, the radio-detected LLAGNs are closer to the low- 
luminosity extrapolation of FR I radio galaxies than to 
that of 'classical' Seyferts. 

Since the Ha luminosity is often used to represent 
the emission-line luminosity for AGNs we now look at 
the dependence of this quantity on the galaxy luminos- 
ity (Fig. ^3|). The Ha-t-[N II] luminosity is significantly 
correlated with galaxy total luminosity separately for all 
ellipticals and all non-ellipticals in the Palomar sample 
(Fig. |l3|b; Table ^). Figure [l^ a suggests that, at a given 
galaxy bulge magnitude, the radio detected ellipticals (ex- 
cept for NGC 4278) have lower emission-line luminosi- 
ties than the radio detected non-ellipticals, though the 
number of ellipticals is only 10. However, when the to- 
tal galaxy luminosity is considered (Fig. |l3|b), the ellip- 
ticals and non-ellipticals detected in the radio fall more 
or less together, i.e. the galaxy total luminosity, rather 
than the galaxy bulge luminosity, appears to be a better 
predictor of the nuclear emission-line luminosity among 
radio detected LLAGNs. The low numerical significance 
of the statistical tests involving the radio-detected galax- 
ies (Table ^) may simply reflect the small number of nuclei 
considered; the correlation between emission-line luminos- 
ity and galaxy bulge or total luminosity becomes signif- 
icant when we also consider the 2 cm non-detections at 
D< 19 Mpc (Table H, bottom). A hnear fit to aU radio- 



The radio power of the radio-detected nuclei which are not 
in ellipticals does not depend on the galaxy total, bulge, 
or disk luminosity; linear fits in all cases give a line with 
approximately zero slope (Fig. |l^ . The upper limits to 
the radio emission of the undetected (small open symbols 
in Fig.|lJ) non-ellipticals strongly support this result. The 
few elliptical nuclei in the sample show a variation of 1 to 
3 orders of magnitude in radio power for a given galaxy 
luminosity. The ellipticals in the sample generally have 
higher radio powers and higher galaxy luminosities than 
the non-ellipticals in the sample. Probably as a result of 
this, when all radio detected nuclei in the sample (i.e. both 
ellipticals and non-ellipticals) are considered together, the 
radio power is found to be correlated with the galaxy bulge 
luminosity (but not the galaxy total luminosity; Table 
with linear fit P2cm ^ Lb ^ (bulge). The slope of this de- 
pendence is much shallower than that seen over 6 orders 
of magnitude of core radio power by Nelson fc Whittle 
(1996) for both powerful radio (FR I and FR II) galaxies 
and Seyfert (type 1 and type 2) galaxies (dotted line in 
Fig. |l|a; P|°J^ oc L|2 (bulge)). It is also shallower than 
the relationship displayed by "normal" nearby elliptical 
and SO gala xiesfl (the 30th per centile of the radio power 
rises as Lg^; Sadler et al., 198E). Figure shows that the 
core radio powers of LLAGNs are less than what would be 
expected from the low-luminosity extrapolations of more 
powerful AGNs. Factors which could affect the slopes and 
radio powers include the differences in the frequencies and 
resolutions between the various datasets. Given the longer 
wavelength (2 cm) and lower resolutio n (1^' — 3") of the 
data u sed by Nelson fc Whittle (1996| ) and Sadler et al 



1989|) , the nuclear radio emission is more likely to include 



|13| gives L 



Ha + [NII]OC l^B 



0.43 



(total) 



detected nuclei in Fig. 
and cx L|^- 26 (bulge). 

As with the [O I] luminosity, the Hq!+ [N II] lumi- 
nosities of the LLAGNs fall closer to the low-luminosity 



emission from arcsec-scale radio jets and star-formation 
activity in the bulge in their data (see e.g. Sect. 11.2 of 
Whittle, 1992b| ) as compared to our higher resolution 2 cm 
data. High resolution radio observations, capable of isolat- 
ing the flat spectrum pc-scale cores in the more powerful 
AGNs, are needed for comparison with our data. 



extrapolation of FR I radio galaxies (Zirbel fc Baum 
1995|f}han to that of 'classical' Seyferts (Fig. |lj). The 
Ha-f [N II] luminosities of the 'classical' Seyferts were de- 
rived by converting the [OIII]A5007 luminosities listed 
in Whittle (1992a) using standard emission- line ratios 
in Seyferts: Ha/H/? =3.1, [O III] A5007/H/3 = 10, and 
[N II] A6583 /Ha = 1 (|Vcillcux fc Osterbrock, l"987|) . As for 



We used the relationships derived by Nelson & Whittle 



(1996| ) and [Sadler et al. (1989| ) after changing Hq to the value 



used in this paper and transforming their 20 cm luminosities 
to 2 cm assuming optically-thin synchrotron emission (spectral 
index —0.7). Both these transformations change the intercept 
of the line but not the slope. 
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8.3. The Interplay Between Core Radio Power and 
Emission-Line Luminosity 

The interplay between radio and emission-line luminosity 
can be better understood in the context of the following 
two factors (discussed above): (1) for all LLAGNs consid- 
ered t ogether, the 2 cm core power fP§°I^ cx: L^^ (bulge)) 



9. Radio Luminosity and Black Hole Mass 

We now investigate the relationship between the nuclear 
radio luminosity and the mass of the nuclear black hole 
in nearby galaxies. When a direct measurement of the 
black hole mass from high spatial-resolution kinematic 



minosity (LHa-|-[Nii]C< Lg^^(bulge)) with increasing bulge 
lumin ositv. as is also the case for "normal" elliptical and 



measurements (Richstone et al., 1998; Gebhardt et al 



increa des more rapidly than the nuclear emission-line lu- 2000|) is not available, we use the relationship found by 



SO gali lxies (^adlcr et al.. 19^89): (2) with increasing disk 
luminosity, the emission-line luminosity increases, but the 
radio power does not. 

A plot of core radio power against 
Ha-H[N II] AA6548,6583 luminosity for all LLAGNs 
at D < 19 Mpc is shown in Fig. |l^. This type of plot 



Merritt fc Ferrarese (2001 ) between the black hole mass 
and the central stellar velocity dispersion, 0-^ . Merritt fc 
Ferrarese (2001) and Gebhardt et al. (2000, who used 



(Te , the luminosity- weighted linc-of-sight dispersion within 
the half-light radius) have shown that this relationship is 
very tight with its accuracy currently limited only by the 
measurement errors in Uc and Mmdo- Values of CTc have 
been taken, in increasing order of preference, from the 



has hfcn — t«ed — te — separate 'radio-loud' — AGNs — ftt 



radio- quiet' ones (e.g. [Xu et al., 1999 ). Also plotted 



the figure are the low-luminosity extrapolations of the 
same relationship for FR I and FR II radio galaxies from 



Zirbel fc Baum (1995 , dashed lines), a nd the sample of 



'classical' Seyferts from Whittle (1992a , small triangles). 
The values quoted in these papers have been converted 
to a Hubble constant of 75 km s^^ Mpc^^, and the 
Seyfert Hq!-)-[N II] luminosities have been derived from 



the [OIIIJA5007 luminosity as described in Sect. 8.1 



Though there is considerable scatter, factor (1) above 
causes a general trend for high bulge-mass LLAGNs to 
occupy the top center of the plot and low bulge-mass 
LLAGNs to occupy the lower left. However, factor (2) 
allows many of the galaxies with low bulge luminosities 
to sit to the right of their compatriots with higher bulge 
luminosities by dint of their higher disk luminosities. 

It is very interesting that most of the radio-detected 
LLAGNs fall in the region of Fig. 
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where the low- 
luminosity extrapolations of the radio-loud (FR I and 
FR II radio-galaxies) a nd radio-quiet (Seyfert galaxies 
from the compilation of Whittle, 1992a ) nuclei intersect. 
There is thus no evidence for a radio-loud - radio-quiet di- 
vide in the LLAGNs. Most of the radio detected LINERs 
and transition nuclei are found close to the low-luminosity 
extrapolations of the FR I and FR II galaxies. Some of 
the LLAGNs are offset to higher emission-line luminosi- 
ties from the FR I and FR II relationships. Essentially 
all of these LLAGNs are in disk-type galaxies, with disk 
gas probably contributing to the emission-line luminos- 
ity. The operational distinction between 'classical' Seyferts 
and low-luminosity Seyferts is log [(Lhq+[nii])/W] ~ 33, so 
it is not surprising that low-luminosity Seyferts near this 
luminosity fall near the low luminosity extrapolation of 
'classical' Seyferts. It is notable in Fig. [ij that the ellipti- 
cal galaxies tend to have a higher radio power at a given 
emission-line luminosity than the SO or disk galaxies, as 
we already inferred from Fig. ^ 



Lyon-Meudon Extragalactic Database (LEDA; Paturel et 
al., 1997), Prugniel et al. (1998) and its update in the 
Hypercat' database^ , Ferrarese & Merritt (2000) and 



Merritt fc Ferrarese (2001). Ferrarese fc Merritt (2000) 



and Merritt fc Ferrarese (2001) list the value of Cc after 
correction for aperture size. They find that the aperture 
corrections for the galaxies in their study are of order a few 
km s~^ except for some more distant (^60 Mpc) galax- 
ies. Therefore, since most of the galaxies we consider be- 
low are relatively nearby (typically < 30 Mpc), we have 
not made any aperture corrections to the values of CTc 
taken from other sources. In support of this, a compar- 
ison between the values of CTc from Ferrarese fc Merritt 
and those taken from LEDA or Hypercat agree to within 
20% (Fig. ^6|). There is no evidence that the difference in 
the values of CTc increases with galaxy distance or for later 
morphological types. However, the difference between the 
two estimates may increase with decreasing galaxy bulge 
luminosity (Fig. |l6|c). 

In order to minimize confusion by radio emission from 
larger scale regions, e.g. jets/lobes or star formation in 
the bulge, we use mainly sub-arcsecond scale radio fluxes, 
and treat all radio flux measurements made at resolution 
> 150 mas as upper limits to the 'core' radio flux. One 
may expect that even the < 150 mas-resolution data (typ- 
ically < 15 pc) overestimates the flux from the central 
few hundred Schwarzschild radii. However, free- free and 
synchrotron self-absorption may be significant on sub-pc- 
scales, so the central flux could be significantly attenuated 
even at 2 cm (e.g. Ulvestad et al., 1999). Most of the mea- 
surements are made at frequencies of 15 GHz or higher 
and thus avoid some of the problems of contamination by 
radio outflows and star formation related processes. For 
the handful of measurements at frequencies lower than 
15 GHz, we converted to estimated 15 GHz fluxes using a 
spectral index of 0, which is the median spectral index 
observed for mas-scale core radio emission in LLAGNs 



( Nagar et al., 2001 ). The lowest observed frequency we 
use is 1.4 GHz (20 cm); thus using a spectral index in the 



littp: / /www-obs.univ-lyonl .fr/hypercat / 
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range ±0.7 instead of would change the extrapolated 
15 GHz flux by a factor of <5. 



9.1. Results for the Distance-Limited Sample of 
LLAGNs 

The correlation between radio luminosity and bulge lumi- 
nosity for LLAGNs has been explored in Sect. The 
basic result is that the two quantities are correlated if all 
radio detections are considered together, but not when 
only non-elliptical galaxies are considered. The plot of 
radio luminosity vs. MDO mass (Fig. ^7|) also shows a 
strong correlation but considerable scatter. In fact the five 
massive ellipticals in the figure - NGC 4278, NGG 4374, 
NGC 4472, NGG 4486, and NGG 4552 - have very simi- 
lar MDO masses (~ 10^ M©) but widely different radio 
powers, even at mas resolution. This large spread does not 
necessarily imply a large spread in the radio luminosities 
within a few hundred as there is evidence that the ra- 
dio emission in all five of these sources is jet dominated 



(see Nagar et al., 2001) 



9.2. Results for a Larger Sample of Nearby Galaxies 

We now consider the interplay between radio power, MDO 
mass and bulge luminosity in the larger sample (Sect. 2.2) 
of relatively nearby galaxies. The core radio fluxes for the 
Seyferts and LINERs at D > 19 Mpc in the Palomar sam- 
ple were taken from Nagar et al. (2000| ) and Nagar et al. 
(in preparation). These are 150 mas resolution, 2 cm ob- 
servations, like the observations of the D < 19 Mpc sample 
of LLAGNs. The core radio fluxes of the additional galax- 
ies were take n from [Garral, Turner, fc Ho (19901 ) iGondon l 
et al. (1990|), pranc (19791), [Crane et al. (19931), |Dhaw an 



KeUerman, fc Romney (1998|), |Di Matteo et al. (2001|) 



Fabbiano, Gioia fc Trinchieri (1989| ), 



1987), Franceschini et al. (1998), Gallimore, Baum fc 



Feretti fc Giovannini 



0'Dea | (1997D , jGcldzalilLT fc FumaluuL (19841), i Gluvamilm 



et al. ( 199^, [ Gralmm el al. (198l|), [ Juucs (1974|), |juucs 



fc Wehrle (19971), [Kellermann et al. (1998|), Krichbaum 



et al. (19981), [Laurent-Muehleisen et al. (19971), jNagar et 



al. (1999|), ^adler et al. (1995|), |Slee et al. (1994|), [Trotter 



and Wrobel fc Heeschen (1991 ). When kinematically- 



et al (199^), [Venturi et al. (1993[ ), [White et al. (1997[ ) 



determined MDO masses were not available, we estimated 
the MDO mass from ctc, as outlined in the previous sec- 
tions. When available, distances have been taken from 
the same reference as that providing the black hole mass. 
Galaxy morphological types, corrected total magnitudes 
in the B band, and distances for the remaining galaxies 
have been taken from LEDA. All quantities were converted 
to the value of Hq - 75 km s^^ Mpc~^ - used in this paper. 
In all, we found that 154 of these galaxies, 37 of which are 
ellipticals, have values for all three of the following: esti- 
mated MDO mass, core radio flux, and bulge luminosity. 
Sixty six of these 154 galaxies have a radio core detected 
at resolution better than 150 mas. 



The nearby galaxies, spanning 9 orders of magnitude 
in radio power, 6 orders of magnitude in MDO mass, and 
5 or ders o f magnitude in bulge magnitude, are plotted in 
Fig. 1S^°. With such a large range of values, it is clear 



that all three quantities are correlated. Asurv tests on 
all 154 galaxies and on only the 66 with radio cores at 
< 150 mas resolution, indicate that the radio power is 
correlated with both the MDO mass and the bulge lu- 
minosity at the 99.99% signiflcance level (Table ||). For 
our Galaxy, we plot the radio flux both at 10 mas reso- 
lution and at ~ 10' (23 pc) resolution; the latter is well 
matched to the linear resolution for the other sources. At 
this matched resolution (which includes emission from the 
whole Sgr A complex), the Galaxy is no longer an outlier 
in the correlations. 

Visually, the relationship between radio power and 
bulge luminosity appears to have less scatter than 
that between radio power and MDO mass, mainly be- 
cause of the outliers in the latter: the massive ellipti- 
cals at the highest radio powers and NGG 3115. For 
Mmdo — 10*'^^^, there is a range of a factor of lO'^ in 
P2cm- We compare the correlations between the three 
variables using the Kendall t coefficient. Numerically, 

-PconcordanceZ-Pdiscordancc = (l + r)/(l-T), whcrC P is the 

probability. That is, r = 1 implies perfect concordance 
(i.e. perfect proportionality), r = — 1 implies complete dis- 
cordance (i.e. perfect inverse proportionality), and t — 
implies an absence of correlation. If we consider all the 
galaxies in Fig. |l^, the Kendall r coefficient indicates 
that the closest correlation is between black hole mass 
and bulge luminosity (Kendall's t— 0.53), while the radio 
power is more or less equally correlated with the black 
hole mass (t= 0.32) and bulge luminosity (r= 0.29). The 
significance of all three correlations is 99.9% (Table |5|). In 
order to determine the primary correlation between the 
three quantities (radio luminosity, black hole mass and 
bulge luminosity), we use a partial correlation analysis 



based on the Kendall r coefficient (Akritas fc Siebert 



1996). This method determines the partial correlation be- 



tween the independent and dependent variables after tak- 
ing away the influence of any correlation with a third (test) 
variable. We get similar results for both cases - partially 
correlating Mmdo and radio power with L (bulge) as the 
test variable (partial t— 0.20, signiflcance 95%), and par- 
tially correlating L(bulge) and radio power with Mmdo as 
the test variable (partial r= 0.16, signiflcance 95%). The 
presence of a 95% signiflcance for the partial correlations 
indicates that the radio power is dependent on both the 
MDO mass and the bulge luminosity. This, however, is 
only a ~2cr result. Thus, at this point, we cannot deter- 
mine whether the radio power is primarily related to the 
black hole mass or the bulge luminosity or to both equally. 
A more conclusive analysis must await a better determi- 



NGC 4278 has moved significantly betwe en Figs. [l7| and 
as H97a use a distance of 9.7 Mpc while Richstone et al 



(1998[ ) use a distance of 17.5 Mpc, which is used in Fig. [l8|. 



14 



Nagar, N. M. et al.: AGNs in a Distance-Limited Sample of LLAGNs 



nation of the MDO versus CTc relationship in lower mass 
ellipticals and in non-ellipticals. 

Among the 37 ellipticals the radio power is correlated 
with the bulge luminosity at the 98.5%-99.5% significance 
level and with the MDO mass at the 98.5%-99.8% level 
(Table 1^). Linear regression analyses by the Buckley James 
method in Asurv yields: 

log Pacnf '''' = 2.33(±0.67) log Lb (bulge) - 4.35 

log Ptl^'''"' = 2.11(±0.63) log Mmdo + 2.33 

Here, and in the equations below, P2cm is measured in 

Watt Hz~^, Lb in solar luminosities, and Mmdo in solar 

masses. 

For the 48 non-elliptical galaxies which are detected 
in the radio at resolutions better than 150 mas, the radio 
power is correlated with both the MDO mass and the 
bulge luminosity at the 99.99% with the following linear 
regression fits 

log Pa™"" = L92(±0.27) log LB(bulge) + 1.14 

log PtlT"" = L15(±0.21) log Mmdo + 11.35 

The presence of a few upper limits to the MDO mass 
and/or bulge luminosity among the non-elliptical galaxies 
which are not detected in the radio complicates the 
statistical tests. When both variables have upper limits 
to some data points, only the Schmitt binned regression 
can be used for linear fitting within the Asurv package. 
The results in this case are quite dependent on the 
number of bins chosen, e.g. for all galaxies with both 
MDO estimates and radio powers, P2cm « M^|jq when 
5 bins are used, and M^^q when 6 bins are used. We 
therefore delete the four nuclei with MDO mass upper 
limits and the three nuclei with upper limits to their 
bulge luminosities from the sample. Having done this, 
the Buckley James linear regression method gives the 
following for all (i.e. whether or not they are detected in 
the radio) non-elliptical galaxies: 
log p^ii^so/iatcr ^ ^ 56(±o.30) log Lg (bulgc) + 2.85 
log Pf^^"^'''"'' = 1.07(±0.21) log Mmdo + 10.52 



The standard deviations of the two relationships are 
similar for the ellipticals and non-ellipticals. Even though 
the linear regression gives a smaller slope for the non- 
ellipticals than the ellipticals in both relationships, this 
result is not certain given the large number of radio upper 
limits at the lower MDO masses and bulge luminosities. 
These data points could steepen the relationship for non- 
ellipticals and bring it closer to that for ellipticals. In fact, 
the current data show no contradiction to all of the galax- 
ies being described by the same relationship. If we do this 
- i.e. fit ellipticals and non-ellipticals together - we get: 
log PfL = 1.89(±0.21) log LB(bulge) - 0.17 
log PfL = 1.31(±0.16) log Mmdo + 8.77 
if all galaxies are considered; and 

1.72(±0.20) log LB(bulge) + 3.15 



li^rr TDradio core 
iOg i^2cm 



log P^2cm° = 1.14(±0.16) log Mmdo + 11-49 
if only those galaxies with radio cores detected at resolu- 
tion < 150 mas are considered. 



10. Discussion 



For all 16 LLA GN s at D < 19 Mpc observed at mas res- 
olution (Sects. L3, 5A) four factors point strongly to a 
non-thermal origin of their mas-scale radio cores: (1) the 
brightness temperatures at 6 cm are too high to be ex- 
plained by thermal processes; (2) the morphologies of the 
five LLAGNs with extended mas-scale radio structure are 
suggestive of pc-scale "jets"; (3) the spectral shapes of 
the radio cores support a non-thermal origin for the ra- 
dio emission (Nagar et al., 2001, 2002); and (4) significant 
flux variability is observed and is much more likely in non- 
thermal than thermal sources. Thus, the mas-scale radio 
emission is almost certainly related to accretion onto a 
massive black hole. The high-brightness temperatures of 
the cores and the presence of either a single component, 
or jet structures about a single bright component, make it 
likely that the radio cores we detect are within ~ 1 pc of 
the accreting black hole. If we are indeed looking directly 
at the accreting region or the base of the jet in the radio, 
then these nuclei are unique among nearby galaxies and 
merit further study in the radio. 

The radio results imply that a large fraction (perhaps 
all) of LLAGNs have accreting massive black holes. If we 
consider only the mas detections, then at least 25% ± 6% 
of LINERs and low-luminosity Seyferts have accreting 
black holes. VLA-detected compact radio cores with flux 
< 2.7 mJy were not investigated with the VLBA; in other 
respects these cores are similar to those with detected mas 
scale structure. Thus it is likely that all LLAGNs with 
VLA-detected compact radio cores (38% ± 8% of LINERs 
and low-luminosity Seyferts) have accreting black holes. 
The scalings between radio power, emission-line luminos- 
ity, and galaxy luminosity provide evidence that the radio 
non-detections are simply lower power versions of the ra- 
dio detections. In fact we find no reason to disbelieve that 
all LLAGNs have an accreting black hole. 

One of the most important results in this paper is that 
compact radio cores are found almost exclusively in mas- 
sive ellipticals and type 1 LLAGNs. For massive ellipti- 
cals, the high bulge luminosity and black hole mass ap- 
pear to be key factors related to the production of a radio 
core, in light of the scalings seen between radio power and 
these parameters. Among non-ellipticals, the preferential 
detection of type 1 LLAGNs may result from the limited 
sensitivity of optical and radio observations, which de- 
tect broad Ha and radio cores in only the more luminous 
LLAGNs. For example, it may be that type 1 LLAGNs 
are in an outburst phase, during which accretion power 
contributes significantly to the total nuclear power out- 
put, and during which they temporarily host both broad 
Ha emission and a compact radio core. This scenario is 
supported by the findings that the observed 2-10 keV lu- 
minosity in type 2 LLAGNs is insufficient to power the op- 
tical emission-lines unless the X-rays are heavily absorbed 
( Tcrashima, Ho, fc Ptak, 2000D , and that the nuclear UV 
emission in some type 2 LINERs, the far-UV component 
of which is responsible for powering the emission lines, is 
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dominated by massive-star clusters ( Maoz ct al., 199^ . 
As another alternative, one can invoke the unified scheme 



( Antonucci, 1993) and posit that all LLAGNs have accret- 
ing black holes and either (a) the radio emission in type 1 
LLAGNs is beamed (weakly relativistic jets [7 ^ 2] can 
give boost factors of up to ~ 5) and/or (b) the 2 cm ra- 
dio emission in type 2 LLAGNs is free-free absorbed by a 
'torus' i.e. tisqhz > 1- 

For the distance-limited (D < 19 Mpc) sample of 
LLAGNs, we found (Sect. 8.3) that the emission-line 
luminosity increases with the galaxy disk luminosity, 
while the radio power is primarily correlated with the 
galaxy bulge luminosity. Thus doubling the bulge luminos- 
ity almost doubles the Prad/Lcmission-iinc ratio (Sect. H, 
Figs. [T^) The lower bulge luminosity nf late-t^^pe 



In the context of a jet model, the correlations between 
core radio power and both black hole mass and galaxy 
bulge luminosity imply either (a) the strength of the pc- 
scale radio jet scales with the mass of the black hole, 
and/or (b) the radio emission from the pc-scale jet is in- 
fluenced by the potential of the galaxy bulge. In the case 
of (a), if the radio power is linearly dependent on the black 
hole mass, as is consistent with Fig. 18, then (following the 
arguments of Ho, 200^ ) one can estimate the radio lumi- 
nosity empirically from the bolometric luminosity (Lboi) 
of the AGN. This method implies that the nuclei with 
detected radio cores in Fig. 



18 have Leoi/LEdd 



10- 



igs. 

galaxids impUpt; hath a smallor raHin pm»rpr nnrl 



to 10~^, with the linear fit to the radio detections (solid 
line in upper panel of Fig. |l8|) lying approxi mately at 
LRn i/LEdd ~ 10 ~'^. The maximal jet model of Falcke fc 

Ti^-^: to fit the 



(19961) predicts that P^°^<; 



smaller Piad/Lcmission-iinc ratio than found for ellipti- 
cals. Thus with increasing total galaxy luminosity, the 
Prad/Lcmission-iinc ratio of elliptical and non-elliptical 
galaxy nuclei will become increasingly separated and can 
create a trend not unlike the 'radio-loud' and 'r adio-quiet 



sequences seen in more powerful AGNs (e.g. |Xu et al 
19990 

In our larger, 'MDO' sample (Sect. |2.2|), the radio 



power is correlated with both the estimated MDO mass 
-"^^ cx M^no) ^'^d the luminosity of the galaxy bulge 



(P 



all 
2cm 



Lg (bulge)) over five or more orders of magni- 



tude (Fig. [18). At this point it is not clear which of the 



Biermann 

slope of the data in Fig. |l^ this requires m to scale roughly 
linearly with Mmdo , equivalent to requiring the LLAGNs 
to accrete at a roughly fixed fraction of the Eddington 
rate. If m does scale with Mmdo this model can reproduce 
the observed radio luminosities of the nuclei in Fig. |l8|even 
for LBoi/LEdd 10"^ (comparing Fig. |l^ with Fig. 1 of 
Falcke & Biermann and assuming their Ldisk ^ Lboi)- 

In the context of (b) above, the relationship between 
radio power and black hole mass may come about indi- 
rectly as outlined below. The core radio power, at arcsec- 
scale resolution, is strongly correlated with the bulge mag- 
nitude over 6 orders of radio power for FR I and FR II 



two is i thc primary correlation. Partial correlation analy- 



radio galaxies, and Seyfert galaxies (Nelson & Whittle 



sis sug gests (at the ' ^ 2i7 level) that both correlations arc 
physically significant. If so, the latter correlation could be 
due to the influence of the bulge mass on the accretion 
rate. More complete results must await a better charac- 
terizatiuii of Lhu 
elliptic ils and iii liuii-ulllpLicals 



1996). For several reasons, powerful radio galaxies and lu- 



minous Seyfert galaxies are not posited to be powered by 
low-radiative-luminosity models, e.g. their derived accre- 
tion rates can be m ^ 10~^'^, which is the currently be- 



f LliL! MDO vcisiis Uc lelaLiuuship lu luw mass 



lieved upper limit for an ADAFs existence ( Narayan et al 



One model which may account for the correlation 
between radio power and black hole mass posits that 
the radio emission is thermal synchrotron radiation from 
the accretion inflow in an ADAF or CDAF type model. 
The approxi mate proportionality between these two quan- 
tities (Sect. |9.2|) is consistent with the L2cm 3 x 
ergs relationship expected in the 



1998), an d they show rapid X-ray variability (e.g. Ptak ct 
al., 1998). Thus, it is possible that some other factor is 



1036 A ff ' mgf 



responsible for the scaling of the core radio power with 
the host galaxy bulge luminosity. This factor could be the 
mass loss rate of stars in the bulge, which should be re- 
lated to both the black hole accretion rate and the pressure 
of the interstellar medium in the vicinity of the nucleus. 
A higher thermal pressure around the radio sourc e would 
tend to confine it and enhance the radio power (cf. Nelson 



ADAF model of Yi fc Boughn (1999) for values of m fc Whittle, 1996 ). Alternatively, the correlation with black 



10^ to ~ 10^^ of the Eddington accretion rate. Here 
m_3 = m/10~'^, m is the mass accretion rate in units 
of the Eddington rate, and M7 is Mmdo/10''. A strong 
argument against all the radio emission coming from an 
ADAF or CDAF is that many of the galaxies considered 



hole mass could be the primary one. As we have shown, 
the nearby galaxies we consider in Fig. ^ (lower panel) 
follow P2°m oc Lg^ (bulge). Also the black hole mass scales 
with the bulge mass as Mmdo = 0.0 05 Mbuigc, wher e 
Mbuigc = 2.5 X 10^ (Lbuige/10^Lo)i-2 dRichstone eTal 



here sl ow morphological evidence for pc-scale jets, have 
core radio spectra which are flatter than what is expected 



1998). Combining these equations gives P^! 



«m1^o. 

The slope of this relationship is exactly the same as that 



core 
cm 



in simple ADAF or CDAF models ( |Nagar et al., 2001] ) 
and/or turn over in the 10-30 GHz range (Di Matteo et al. 
2001; Nagar et al., in preparation, see Nagar et al. 2002). If 
a significant (perhaps dominant) component of the radio 
emission is from jets, then the data in Fig. ^ imply accre- 
tion rates much lower than ^ 10~^ — 10"*^ of the Eddington 
rate within an ADAF or CDAF type model. 



of the ADAF model of |Yi fc Boughn (1999|) and is not too 
far from the behavior of the data in Fig. |18| (upper panel) . 

11. Conclusions 

The completion of our radio survey of the 96 nearest (D < 
19 Mpc) LLAGNs from the Palomar sample of nearby 
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bright galaxies has yielded the following main results: 

a) All (16 of 16) of the LLAGNs investigated at mas res- 
olution with the VLBA have pc-scale cores with bright- 
ness temperatures 10® K. The five of these nuclei with 
the highest core flux all have pc-scale jets. The luminos- 
ity, brightness temperature, morphology, and variability 
of the radio emission all argue against an origin in star- 
formation related processes or as thermal emission. Thus, 
the core radio emission probably originates either in an ac- 
cretion inflow onto a supermassive black hole or from jets 
launched by this black hole-accretion disk system. The 
latter explanation is supported by the morphologies and 
spectra of the radio cores; 

b) there is no reason to believe that the remain- 
ing LLAGNs with compact radio cores (investigated at 
150 mas resolution) are different from the 16 LLAGNs 
investigated at mas resolution. Thus, at least half of all 
LINERs and low-luminosity Seyferts probably contain ac- 
creting black holes. The incidence for transition nuclei is 
much lower; 

c) compact radio cores are preferentially found in massive 
ellipticals and in type 1 nuclei, i.e. nuclei in which broad 
Ha emission is present. For these nuclei, the core radio 
power is proportional to the broad Ha luminosity. The 
preferential detection of type 1 nuclei could result from: 
1) observational selection effects, in which broad Ha and 
radio cores have been found only in the more powerful 
LLAGNs, 2) only the type 1 LLAGNs are bona-fide AGNs, 
or 3) within the unified scheme, type 1 nuclei are beamed 
and/or type 2 nuclei are free-free absorbed in the radio; 

d) the radio power of the compact core is correlated with 
both the galaxy luminosity and the luminosity and width 
of the nuclear emission lines. These trends suggest that we 
have detected only the brighter LLAGNs, i.e. the true in- 
cidence of accreting black holes in LLAGNs is higher than 
indicated by our radio detections alone. 

e) The core radio and nuclear emission-line properties of 
LLAGNs fall close to the low-luminosity extrapolations 
of more powerful AGNs, providing further support for a 
common central engine; 

f) low-luminosity Seyferts and LINERs share many of 
the same characteristics in the radio. The transition nu- 
clei detected are those which are the closest, in terms of 
emission-line diagnostic ratios, to Seyferts and LINERs. 
Thus at least some transition nuclei are really composite 
Seyfert/LINER + HII region nuclei, with the core radio 
power dependent on the Seyfert/LINER component; 

g) with increasing bulge luminosity, the radio power in- 
creases more rapidly than the emission-line luminosity. 
Also, there is evidence that the disk luminosity of the 
galaxy is correlated with the nuclear emission-line lumi- 
nosity (but not the core radio power). Both these factors 
contribute to a radio loud/quiet division between ellipti- 
cals and non-ellipticals at high bulge luminosities; 

h) investigation of a sample of ~ 150 nearby bright 
galaxies, most of them LLAGNs, shows that the nuclear 
(<150 mas) radio power is strongly correlated with both 
the estimated MDO mass and the galaxy bulge luminosity. 



Partial correlation analysis on the two correlations yields 
the result that each correlation is meaningful even after 
removing the effect of the other correlation. The nature of 
these important correlations is discussed. Low accretion 
rates (< 10~^-10^'^ of the Eddington rate) are implied in 
both ADAF- and jet-type models; 

i) about half of all LLAGNs investigated show significant 
inter-year variability at 2 cm and 3.6 cm. 
In short, all evidence points towards the presence of ac- 
creting black holes in a large fraction, perhaps all, of 
LLAGNs. 
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Fig. 1. Histograms of the galaxy distance for the (a) distance-limited (D< 19 Mpc) LLAGN sample, and (b) MDO 
sample. Virgo cluster members create the peak at D-^ 17 Mpc. 
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Fig. 2. The inter-year variability of the nuclear (left) 2 cm flux, and (right) 3.6 cm flux, in LLAGNs. The Y axis is 
the logarithm of the average measured flux, and the X axis is the r.m.s. of the flux variation about the average flux 
divided by the average flux. For the 3.6 cm data, and a few 2 cm datapoints, only two epochs are available and for 
these cases AS represents the difference of the two flux values. Estimated ±2 a errors in X and Y are shown. When 
error bars in Y are not shown, they are smaller than the plot symbol. Points outside the shaded area are likely to be 
dominated by errors in flux-calibration and/or radio imaging (see text); these points are therefore not reliable. 
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Fig. 3. 2 cm light curves for five of the LLAGNs in Fig. The 2 cm flux at the last epoch (December 2000) has 
been obtained by scaling the previous epoch's 2 cm flux by the ratio of the 3.6 cm fluxes, which were obtained on 
both epochs. Estimated ±2tT error bars (from flux calibration and mapping errors) are shown for each datapoint. In 
each panel the horizontal solid line shows the average flux over the observed epochs and the dashed lines show the 
estimated ± 2 cr errors for the average flux value. 
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Fig. 4. 5 GHz (6 cm) VLBA maps of NGC 4374 (left) and NGC 4552 (right). The contours are integer powers 
of v^, multiplied by the ~ 3cr noise level of 1.2 mJy for NGC 4374, and by the ^ 2a noise level of 0.8 mJy for 
NGC 4552. The peak flux-densities are 152.6 mJy/beam and 93.8 mJy/beam, respectively. 
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Fig. 5. Relationship between [O I] A6300A luminosity and 2 cm core (150 mas resolution) radio power for (a) all 
LINERs and low-luminosity Seyferts, and (b) all transition nuclei, in the distance-limited (D < 19 Mpc) sample. Large 
symbols are used for radio detections. Objects with non-photometric [O I] A6300A flux measurements are shown as 
open circles. 
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Fig. 6. Relationship between the FWHM of the nuclear [N II] A6583 emission-line (as listed in H97a) and the 
inclination-corrected rotational velocity of the galaxy {AV^^^; as listed in H97a) for (a) LINERs, (b) Seyferts, and 
(c) transition nuclei, in the distance-limited (D < 19 Mpc) sample. Radio-detected ellipticals are shown with large 
crosses, radio-detected non-ellipticals with small crosses, and all radio non-detections with dots. 
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Fig. 7. Relationship between the ratio of the FWHM of the nuclear [N II] A6583 emission-line to the inclination 

corrected rotational velocity of the galaxy and the 2 cm core (150 mas resolution) radio power for (a) radio detected 
nuclei and (b) radio non-detected nuclei in the distance-limited (D < 19 Mpc) sample. Large symbols are used for 
elliptical galaxies and filled symbols for radio detections. The symbols in the right panel thus represent upper limits 
to the radio power. 
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Fig. 8. Histograms of (a) the galaxy morphological type, and (b) the local galaxy density, comparing the D < 19 Mpc 
LLAGNs detected at 2 cm (solid lines) and those not detected at 2 cm (dotted lines). All values are as listed in II97a. 
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Fig. 9. (a) The luminosity of the broad component of Ha as a function of the 2 cm core (150 mas resolution) radio 
power for all (20) type 1.9 AGNs in the Palomar sample for which photometric broad Ha fluxes are available. LINERs 



are plotted as circles and Seyferts as triangles. Nuclei with definite and probable detections of broad Ha ( Ho ct al 
1997c0are shown with filled and open symbols, respectively. Elliptical nuclei are shown with large symbols; (b) same 
as (a) but for fiuxcs instead of luminosities. 
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Fig. 10. The ratio of 2 cm core (150 mas resolution) radio luminosity to [OIJA6300 luminosity as a function of 
galaxy bulge magnitude in the blue band for type 1 LLAGNs (a) and type 2 LLAGNs (b) in the distance-limited 
(D < 19 Mpc) sample. LINERs are plotted as circles, Seyferts as triangles, and transition nuclei as squares. Large 
symbols are used for elliptical galaxies and filled symbols are used for nuclei detected at 2 cm. Open symbols therefore 
represent upper limits to the true radio to emission-line luminosity ratio. Nuclei with upper limits for both radio and 
[O I] luminosities are not plotted. 



24 



Nagar, N. M. et al.: AGNs in a Distance-Limited Sample of LLAGNs 



5a, 



o 



OX) 

o 



1.0 



0.5 



0.0 



-0.5 



Seyfert 



oo o o " • 



• y- 

□o □ o 
□ o 



□LINER+Trans 



-0.4 -0.2 0.0 0.2 0.4 0.6 
log ([N II] X6583 / Ha) 
(a) 



5a. 



o 



LD 



OX) 

o 



1.0 



0.5 



0.0 



-0.5 



Seyfert 



o 

b- 



o • 

o 



□ □-□-■-ao---*-o'*--« 



■ % W<* •*• oo 



□ □ 
□ 



° .0 ff 

o on 



LINER+Trans 



-0.6 -0.4 -0.2 -0.0 0.2 0.4 
log ([S II] X6716 / Ha) 
(b) 



1.0 



5a. 



o 



■< 0.5 



OX) 



0.0 



-0.5 



Seyfert 



o o • 

o o o 



»• • 



□ 



□ :o 

CO 



LINER 



-1.0 -0.5 0.0 

log ([0 I] X6300 / Ha) 
(c) 

Fig. 11. Plots of standard diagnostic emission- line ratios 
for the D < 19 Mpc LLAGNs. For clarity, nuclei with upper 
limits to any one of the four emission-lines in each respec- 
tive panel arc not plotted. Emission-line fluxes have been 
taken from H97a. LINERs and low-luminosity Seyferts 
are plotted as circles and transition nuclei are plotted 
as squares. Filled symbols arc used for nuclei detected at 
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Fig. 12. The dependence of the luminosity of the [O I] A6300 emission-line on the absolute magnitude in the blue band 
of the galaxy bulge for the distance-limited (D < 19 Mpc) sample. LINERs are plotted as circles, Seyferts as triangles 
and transition nuclei as squares. Large symbols are used for elliptical galaxies and filled symbols for radio detected 
nuclei. Also shown is the low-luminosity extrapolation of the linear fit to the same relationship for the luminous 
Seyferts in the Whittle (1992a ) compilation (dotted line; see text) and an indicative low -luminosity extrapolation of 
the linear fit to the same relationship for FR I radio galaxies from Zirbel fc Baum {^99^ , dashed line; see text). 
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Fig. 13. The dependence of the luminosity of the Ha-|-[N II] AA6548,6583 emission-lines on the absolute magnitude 
in the blue band of (a) the galaxy bulge; (b) the entire galaxy; and (c) the galaxy disk, for the distance-limited 
(D < 19 Mpc) sample. Some of the galaxies fall beyond the left and bottom boundaries of the plots. Symbols are as 
in the previous figure. Linear fits to the radio-detected nuclei are shown as the solid line in each panel. Also shown in 
each panel is the low-luminosity extrapolation of a linear fit to the LfjQ,+[Nii] vs Mb (bulge) relationship for luminous 
Seyferts in the Whittle (1992) compilation (dotted line; see text) and FR I radio galaxies from Zirbel & Baum (1995, 
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Fig. 14. As for Fig.|T^, but with the 2 cm core (150 mas resolution) radio power as the y axis. Some of the galaxies 
fall beyond the left and bottom boundaries of the plots. Each panel also shows the best linear fit to the radio detected 
nuclei (solid line), and the low- luminosity extrapolation of a linear fit to the core radio power vs bulge luminosity 
relationship in both luminous Seyferts and (FR I and FR II) radio galaxies (Nelson fc Whittle, 1996, dotted line). 
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Fig. 15. A plot of 2 cm core 
(150 mas resolution) ra- 
dio power versus nuclear 
Ha-h[N II] AA6548,6583 lumi- 
nosity for all LLAGNs in the 
distance-limited (D < 19 Mpc) 
sample. Symbols are as used 
in Fig.0. "Classical" Seyfert 
galaxies (from Whittle 1992a) 
are plotted as the smallest 
triangles (for these the [O III] 
luminosity was converted to 
an Ha+[N II] luminosity using 
standard flux ratios for Seyferts, 
see text). Also shown are the 
low-luminosity extrapolations of 
linear fits to the same relation- 
ship for FR I and FR II radio 
galaxies (dashed lines; Zirbel & 
Baum 1995, see text). 
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Merrit : (2000 ) and Merritt fc Ferrarese (2001 ) is plotted against the (a) morphological type of the galaxy; (b) distance 
to the galaxy; and (c) absolute magnitude in the B band of the galaxy bulge, for galaxies in the Palomar sample. The 
values from these sources agree to better than ±20%. 
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Fig. 17. Dependence of the core 
(150 mas resolution) 2 cm power 
on the mass of the supermassive 
black hole for LLAGNs in the 
distance-limited (D < 19 Mpc) 
sample. For the radio detections, 
large symbols are used for ellip- 
ticals (which are labelled with 
their names) and filled sym- 
bols for the most reliable mea- 
surements of MDO mass (see 
text). All radio non-detections 
are shown by the downward 
pointing triangles. NGC 185 lies 
beyond the bottom left corner of 
the plot. 
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Fig. 18. Upper Panel: The dependence of the core radio power on the mass of the supermassivc black hole for 
nearby galaxies. Only galaxies with all three of radio power, MDO estimate and Lb (bulge) are plotted. Radio flux 
measurements with resolution > 150 mas are plotted as upper limits. Ellipticals are plotted as filled circles and non- 
elliptical radio detections as open circles. For clarity the symbols for non-elliptical nuclei with radio upper limits arc 
not plotted (only the downward pointing arrows are plotted for these). Galaxy names are provided for data points in 
the less crowded areas of the plots. NGC 185, with an estimated MDO mass of 6.8x10^ Mq, lies off the left of the 
plot. The two values for Sgr A* correspond to radio powers within 10 mas and within ^ 10'. In each panel, the solid 
and dashed lines show the linear fits to all (66) nuclei with radio cores detected at < 150 mas resolution, and all (154) 
nuclei plotted, respectively. The linear fits take account of the upper limits in both variables. Lower Panel: same as 
upper panel, but for bulge luminosity in the B band instead of black hole mass. 



Nagar, N. M. et al.: AGNs in a Distance-Limited Sample of LLAGNs 



Table 1 

THE DISTANCE- AND 
RADIO-FLUX-LIMITED LLAGN SAMPLE 



Name 


WLIlcl 


■ ■ 

Activity 


L-'ISL. 






Ncune 


Type 


i,ivipc; 














NGC 2787 




L1.9 


13.0 


a 


NGC3031 


M81 


S1.5 


3.6 


b 


NGC 3718 




L1.9 


17.0 


c 


NGC 4143 




L1.9 


17.0 


c 


NGC 4168 




S1.9 


16.8 


c 


NGC 4203 




L1.9 


9.7 


a 


NGC 4258 


M 106 


S1.9 


6.8 


e 


NGC 4278 




L1.9 


9.7 


a 


NGC 4374 


M84 


L2 


16.8 


c 


NGC 4472 


M49 


S2 


16.8 


c 


NGC 4486 


M87 


L2 


16.8 


d 


NGC 4552 


M89 


T2 


16.8 


c 


NGC 4565 




S1.9 


9.7 


a 


NGC 4579 


M58 


S1.9 


16.8 


a 


NGC 4772 




L1.9 


16.3 


c 


NGC 5866 


M 102 


T2 


15.3 


a 



Note. — This table lists all LLAGNs from the Palo- 
mar spectroscopic sample which have D < 19 Mpc and 
S^tm > 2.7 mJy. 

Note. — Columns are: (1) galaxy name; (2) other 
common name; (3) nuclear activity type as derived by 
H97a. 'L' represents LINER, 'S' represents Seyfert, 
and 'T' represents objects with transitional 'L' + 'H' 
spectra. '2' implies that no broad Ha is detected, '1.9' 
implies that broad Ha is present, but not broad H/3, 
and '1.0', '1.2', or '1.5' implies that both broad Ha 
and broad H/3 are detected, with the specific type de- 
pending on the ratio of the flux in [O III] A5007 to the 
flux in broad -I- narrow H/3 (Osterbrock 1981); (4) dis- 
tance (in Mpc) to galaxy, as Usted in H97a, except for 
NGC 3031, for which we use the distance from Freed- 
man et al. (1994), and NGC 4258, for which we use 
the distance from Herrnstein et al. (1999); (5) source 
for VLBI data, a = Falcke et al. (2000a); b = Bartel et 
al. (1982); c = this work; d = Junor & Biretta (1995); 
e = Herrnstein et al. (1997). 
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Table 2 



NEW 2 CM VLA OBSERVATIONS OF LOW-LUMINOSITY AGNs 



Name 


Activity 


T 


R.A. 


Dec. 


A 


Speak 


Stot 


Dist. 


Log P1°S 


a 


Com. 




Type 




(B1950) 


(B1950) 


(") 


(mJy/ 
beam) 


(mJy) 


(Mpc) 


(WHz"') 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


IC239 


L2:: 


6.0 










<0.9 


16.8 


< 19.48 






IC 1727 


T2/L2 


9.0 










<0.9 


8.2 


< 18.86 






NGC 428 


L2/T2: 


9.0 










<0.9 


14.9 


< 19.38 






NGC 660 


T2/H: 


1.0 










<0.9 


11.8 


< 19.18 






NGC 1055 


T2/L2:: 


3.0 










< 1.8 


12.6 


< 19.53 






NGC 1058 


S2 


5.0 










<0.9 


9.1 


< 18.95 






NGC 2541 


T2/H: 


6.0 










< 1.0 


10.6 


< 19.13 






NGC 2683 


L2/S2 


3.0 










<0.9 


5.7 


< 18.54 






NGC 2685 


S2/T2: 


-1.0 










<0.9 


16.2 


< 19.45 






NGC 2787 


L1.9 


-1.0 


09 14 49.469 


69 24 50.82 


1.6 


11.1 


11.4 


13.0 


20.35 


F 


1 


NGC 2841 


L2 


3.0 


09 18 35.841 


51 11 24.10 


7.4 


1.1 


2.1 


12.0 


19.28 


F 


2,8,9 


NGC 3031 


S1.5 


2.0 


09 5127.313 


69 18 08.15 


6.3 


164.1 


164.8 


3.6 


20.41 


F 


3 


NGC 3368 


L2 


2.0 










< 1.0 


8.1 


< 18.89 






NGC 3379 


L2/T2:: 


-5.0 










< 1.0 


8.1 


< 18.89 






NGC 3486 


S2 


5.0 










< 1.0 


7.4 


< 18.82 






NGC 3489 


T2/S2 


-1.0 










< 1.0 


6.4 


< 18.69 






NGC 3623 


L2: 


1.0 










<0.9 


7.3 


< 18.76 






NGC 3627 


T2/S2 


3.0 


11 17 38.475 


13 15 55.76 


8.3 


1.1 


2.9 


6.6 


18.76 


F 


2,4,8 


NGC 3628 


T2 


3.0 


11 17 40.354 


13 51 46.06 


6.6 


1.5 


2.2 


7.7 


19.03 




7,8 


NGC 3675 


T2 


3.0 










< 1.0 


12.8 


< 19.29 






NGC 3718 


L1.9 


1.0 


11 29 49.910 


53 20 38.47 


2.0 


10.5 


10.8 


17.0 


20.56 


F 


2,9 


NGC 3941 


S2: 


-2.0 










< 1.1 


18.9 


< 19.67 






NGC 3982 


S1.9 


3.0 










< 1.0 


17.0 


< 19.54 






NGC 4013 


T2 


3.0 










< 1.0 


17.0 


< 19.54 






NGC 4051 


S1.2 


4.0 










< 1.0 


17.0 


< 19.54 






NGC 4138 


S1.9 


-1.0 


12 06 58.344 


43 57 48.12 


1.2 


1.5 


1.3 


17.0 


19.71 


F 


2,8 


NGC 4143 


L1.9 


-2.0 


12 07 04.548 


42 48 44.26 


1.6 


9.8 


10.0 


17.0 


20.53 


F 


1,9 


NGC 4150 


T2 


-2.0 










< 1.0 


9.7 


< 19.05 






NGC 4168 


S1.9: 


-5.0 


12 09 44.270 


13 28 59.62 


1.0 


3.0 


3.1 


16.8 


20.01 


F 


2,8,10 


NGC 4203 


L1.9 


-3.0 


12 12 33.943 


33 28 30.48 


1.2 


8.8 


9.0 


9.7 


20.00 


F 


1,10 


NGC 4216 


T2 


3.0 


12 13 21.620 


13 25 38.11 


2.1 


1.2 


1.3 


16.8 


19.61 


F 


2,8,10 


NGC 4258 


S1.9 


4.0 


12 16 29.365 


47 34 53.19 


0.4 


2.6 


3.0 


6.8 


19.16 


F 


2,8,9 


NGC 4278 


L1.9 


-5.0 


12 17 36.159 


29 33 29.27 


1.6 


83.8 


87.7 


9.7 


20.97 


F 


1,10 


NGC 4293 


L2 


0.0 


12 18 40.996 


18 39 34.95 


5.7 


0.7 


1.4 


17.0 


19.38 


S 


6,8,10 


NGC 4314 


L2 


1.0 










< 1.0 


9.7 


< 19.05 






NGC 4321 


T2 


4.0 










<0.9 


16.8 


< 19.48 






NGC 4346 


L2:: 


-2.0 










< 1.0 


17.0 


< 19.54 






NGC 4350 


T2:: 


-2.0 










<0.9 


16.8 


< 19.48 






NGC 4374 


L2 


-5.0 


12 22 31.584 


13 09 49.82 


3.3 


180.7 


183.7 


16.8 


21.79 


F 


2,10 


NGC 4388 


S1.9 


3.0 


12 23 14.641 


12 56 20.09 


2.1 


2.2 


3.7 


16.8 


19.87 




10 


NGC 4394 


L2 


3.0 










<0.9 


16.8 


< 19.48 






NGC 4395 


S1.8 


9.0 










<0.9 


3.6 


< 18.14 






NGC 4414 


T2: 


5.0 










<0.9 


9.7 


< 19.01 






NGC 4419 


T2 


1.0 


12 24 24.716 


15 19 26.53 


2.8 


2.7 


3.6 


16.8 


19.96 


S 


6,10 


NGC 4438 


L1.9 


0.0 










<0.9 


16.8 


< 19.48 






NGC 4450 


L1.9 


2.0 


12 25 58.284 


17 21 40.94 


1.5 


2.0 


2.7 


16.8 


19.83 


F 


2,8,10 


NGC 4457 


L2 


0.0 










< 1.0 


17.4 


< 19.56 






NGC 4459 


T2: 


-1.0 










< 1.0 


16.8 


< 19.53 






NGC 4472 


82:: 


-5.0 


12 27 14.183 


08 16 35.99 


1.9 


3.7 


4.1 


16.8 


20.10 






NGC 4477 


S2 


-2.0 










< 1.0 


16.8 


< 19.53 






NGC 4486 


L2 


^.0 


12 28 17.568 


12 40 01.74 


3.1 


2725.7 


2835.7 


16.8 


22.96 


F 


5 


NGC 4494 


L2:: 


-5.0 










<0.8 


9.7 


< 18.95 






NGC 4501 


S2 


3.0 










< 1.1 


16.8 


< 19.57 






NGC 4548 


L2 


3.0 


12 32 55.275 


14 46 17.67 


2.5 


1.4 


1.6 


16.8 


19.67 


F 


1,8,9 


NGC 4552 


T2: 


-5.0 


12 33 08.293 


12 49 53.58 


2.7 


58.1 


58.6 


16.8 


21.29 


F 


2 
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Table 2 — Continued 



Name 


Activity 


T 


R.A. 


Dec. 


A 


Speak 


Slot 


Dist. 




a 


Com. 




Type 




(B1950) 


(B1950) 


(") 


(mJy/ 

beam) 


(mJy) 


(Mpc) 


(W Hz"') 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


NGC 4565 


S1.9 


3.0 


12 33 52.017 


26 15 45.89 


4.4 


3.1 


3.1 


9.7 


19.54 


F 


1,8,11 


NGC 4579 


S1.9/L1.9 


3.0 


12 35 11.999 


12 05 34.85 


3.1 


20.8 


20.6 


16.8 


20.85 


F 


1 


NGC 4596 


L2:: 


-1.0 










< 1.1 


16.8 


< 19.57 






NGC 4636 


L1.9 


-5.0 


12 40 16.661 


02 57 41.59 


2.6 


1.6 


1.8 


17.0 


19.74 




8 


NGC 4698 


S2 


2.0 










< 1.0 


16.8 


< 19.53 






NGC 4713 


T2 


7.0 










< 1.1 


17.9 


< 19.63 






NGC 4725 


S2: 


2.0 










< 0.9 


12.4 


< 19.22 






NGC 4736 


L2 


2.0 


12 48 31.915 


41 23 31.35 


6.3 


1.9 


1.7 


4.3 


18.62 


F 


2,8 


NGC 4762 


L2: 


-2.0 


12 50 25.139 


11 30 03.27 


2.8 


0.9 


1.3 


16.8 


19.48 


F 


2,8 


NGC 4772 


L1.9 


1.0 


12 50 56.003 


02 26 22.36 


2.5 


3.3 


3.4 


16.3 


20.02 


F 


2,8 


NGC 4826 


T2 


2.0 










< 0.9 


4.1 


< 18.26 






NGC 5194 


S2 


4.0 










< 1.1 


7.7 


< 18.89 






NGC 5195 


L2: 


10.0 










< 1.1 


9.3 


< 19.06 






NGC 5879 


T2/L2 


4.0 










< 1.1 


16.8 


< 19.57 






NGC 6503 


T2/S2: 


6.0 










< 1.0 


6.1 


< 18.65 






NGC 7177 


T2 


3.0 










< 1.1 


18.2 


< 19.64 






NGC 7331 


T2 


3.0 










< 1.1 


14.3 


< 19.43 







Note. — Columns are: (1) galaxy name; (2) nuclear activity type as given by H97a. 'L' represents LINER, 'S' represents Seyfert, 
'H' represents an H II region type spectrum, and 'T' represents objects with transitional 'L' + 'H' spectra. '2' implies that no broad 
Hct is detected, '1.9' implies that broad Ha is present, but not broad H/3, and '1.0' or '1.5' implies that both broad Hq and broad H/3 
are detected, with the specific type depending on the ratio of the flux in [O III] A5007 to the flux in broad + narrow H/3 (Osterbroclc 
I98I). The ':' and '::' symbols represent imcertain, and highly imcertain, classifications, respectively; (3) Hubble morphological 
parameter T as listed in H97a; (4) and (5) 2 cm radio position; (6) offset, in arcsec, between the 2 cm radio position and the position 
of the galaxy optical nucleus listed in Cotton et al. (1999); (7) and (8) peak flux-density and total flux, obtained by fitting a single 
Gaussian (with peak flux-density, major and minor axes as free parameters) to the nuclear radio source in our maps. Fluxes were 
derived from maps with resolution ^ 150 mas for all radio detected sources (see Sect. 4.1). For most sources, the \a error in the flux 
determination comes from a 2.5% flux-bootstrapping error and the 0.2 mJy r.m.s. noise in the final maps. Additional errors in the 
flux determination, if present, are fisted in the conunents column; (9) distance to galaxy, as listed in H97a, except for NGC 3031 for 
which we use the distance from Freedman et al. (1994); (10) the logarithm of the core radio power (derived from the peak radio flux- 
density); (11) indication of whether the radio spectrum of the 2 cm detected core is Steep (a < -0.3; oc v"), or Flat (a > -0.3); 
(12) comments as listed below. 

Note. — Comments; 1; column 1 1 taken from Nagar et al. (2000); 2: column 1 1 derived from comparison with FIRST; 3: column 
1 1 derived from comprison with NVSS; 4: colunm 1 1 derived from comparison with Laurent-Muehleisen et al. (1997); 5: column 1 1 

derived from Nagar et al. (2001); 6: column 11 derived from VLA archive data (see Sect. 4.1); 7: the total 2 cm flux we measured, 
including the extended starburst region, is 28 mJy; 8: accurate self-calibration was not possible, therefore the measured flux may be 
lower than the true flux due to decorrelation losses; 9: additional ~1.5% flux error due to elevation effects. 10: additional ~2.5% 
flux error due to elevation effects. 11: additional ~4% flux error due to elevation effects. 



Table 3 

NEW OBSERVATIONS OF LLAGN WITH THE VLBA 



Name 


Activity 


R.A. 


DEC. 


Peak Flux 


total Flux 


Beam 


Beam 


Beam 


Log Tb 


Com. 




Type 


(J 2000) 


(J 2000) 


(mjy/beam) 


(mjy) 


(mas) 


(pc) 


P.A. 


(K) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 2655 


S2 






<0.76 




4.8x4.3 


0.57x0.51 


45.0 


<6.4 




NGC 3147 


S2 


10 16 53.6509 


73 24 02.695 


6.1 


8.1 


1.6x1.3 


0.32x0.25 


26.8 


>8.3 


a 


NGC 3718 


L1.9 


11 32 34.8530 


53 04 04.518 


3.4 


5.3 


2.2x1.7 


0.18x0.14 


-54.3 


>7.8 


b 


NGC 4143 


L1.9 


12 09 36.0679 


42 32 03.036 


8.1 


8.7 


2.3x1.2 


0.19x0.10 


-41.2 


>8.3 




NGC 4168 


S1.9 


12 12 17.2685 


13 12 18.701 


2.5 


2.4 


3.2x1.6 


0.26x0.13 


10.3 


>7.6 


b 


NGC 4374 (Core) 


L2 


12 25 03.7433 


12 53 13.142 


160.0 


160.0 


2.7x1.4 


0.22x0.11 


6.6 


>9.5 




NGC 4374 (North) 








11.9 


16.6 










c 


NGC 4374 (North2) 








2.4 


2.4 










c,d 


NGC 4374 (South) 








7.0 


7.0 










c 


NGC 4472 


S2 


12 29 46.7619 


08 00 01.713 


2.2 


2.3 


3.4x1.5 


0.28x0.12 


5.2 


>7.5 


b 


NGC 4552 (Core) 


T2 


12 35 39.8073 


12 33 22.829 


99.5 


99.5 


2.9x1.2 


0.23x0.10 


2.4 


>9.3 




NGC 4552 (Eastl) 








7.9 


7.9 










c 


NGC 4552 (East2) 








5.3 


14.6 










c,e 


NGC 4552 (West) 








9.8 


17.0 










c 


NGC 4772 


L1.9 


12 53 29.1613 


02 10 06.157 


1.8 


1.3 


5.0x4.0 


0.39x0.32 





>6.8 


b,f 



Note. — Colurmis are: (1) galaxy name; (2) activity type as derived by H97a and explained in the footnotes to Table 1; (3) and (4) nuclear radio 
position. Positions were determined by first fitting a Gaussian to the central component in the non-self-calibrated VLBA maps, and then correcting the 
measured positions to reflect the more accurate phase caUbrator positions now available (see Sect. 4.3); (5) and (6) peak flux-density and total flux, 
measured, unless noted otherwise, by fitting a Gaussian to the final map; (7) the major and minor axes of the synthesized beam in milUarcsec (mas); 
(8) the major and minor axes of the synthesized beam in parsec using the distance listed in H97a; (9) the position angle (P.A.) of the beam in degrees, 
measured N through E; (10) logarithm of the brightness temperature in Kelvin of the radio core derived using the formula listed in Falcke et al. (2000a). 
Since most of the sources are unresolved, the values from the formula are usually lower limits to the true brightness temperature; (11) comments as listed 
below. 

Note. — Comments: (a) not in the sample listed in Table 1; (b) peak flux-density (from non-self-calibrated map) multiplied by 1.3 to account for 
decorrelation losses (see Sect. 3.2); (c) from a single Gaussian fit; (d) northernmost component in Fig. la; (e) easternmost component in Fig. lb; (f) all 
data from the longest baselines (to antennas on Mauna Kea and Saint Croix), were deleted because of large phase wraps; (g) the peak flux-density listed 
is a 10 (T upper limit, from a naturally weighted map, obtained after deleting all data from the antennas at Mauna Kea and Saint Croix. 
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Table 4 

RESULTS OF BIVARIATE STATISTICAL COMPARISONS FOR THE DISTANCE-LIMITED (D < 19 Mpc) LLAGN SAMPLE 



Sub-sample Variable 


Variable 


Cox 


Kendall 


Spearman 


A 


B 




r 


P 


(1) (2) 


(3) 


(4) 


(5) 


(6) 




All (31) LLAGNs with 2 cm radio cores in the distance-limited (D < 19 Mpc) sample 


All nuclei Mb (bulge) 




98.5 


92.5 


92.5 


Mb (total) 


P2cm 


71.0 


31.0 


21.0 


Mb (bulge) 


Lhc<+[NII] 


73.0 


68.0 


66.0 


Mb (total) 


LHa+[NII] 


92.7 


88.0 


89.0 


Lhc«+[NII] 


P2cm 


98.8 


97.7 


97.0 


L[oi] 


P2cm 


99.99 


99.32 


98.4 


FWHM([N II]) 


L[Nn] 


99.06 


99.47 


99.51 


FWHM([N II]) 


P2cm 


99.99 


99.99 


99.97 


L[NII] 


P2cm 


99.15 


98.5 


97.6 


Mmdo 


P2cm 


99.75 


99.28 


99.10 


P2cm 


FWHM([NII])/AV?ot 


99.97 


99.94 


99.91 


SO and later only (24 galaxies) L[oi] 


P2cm 


99.94 


99.64 


99.30 


LHa+[NII] 


P2cm 


88.81 


93.72 


91.39 


P2cm 


FWHM([NII])/AVrot 


98.3 


99.20 


99.29 


As above, plus all LINERs and Seyferts not detected at 2 cm in the distance-limited (D < 19 Mpc) sample 


All nuclei Mb (bulge) 


P2cm 


99.88 


99.82 


99.87 


Mb (total) 


P2cm 


99.54 


94.9 


98.2 


LHa+[Nn] 


P2cm 




99.16 


99.49 


L|oi 


P2c,„ 




99.97 


99.97 


All LLAGNs not ciclcclcd al 2 cm in Ihc dislancc-liniilcd (1) < 19 Mpc, 


1 sample 




All nuclei Mb (bulge) 


LHa+[NII] 


99.99 


99.99 


99.99 


Mb (total) 


LHa+[NII] 


99.87 


99.99 


99.99 


All LLAGNs in the distance-limited (D < 19 Mpc) sample 


All nuclei Mb (bulge) 


Lhq+[NII] 


99.98 


99.99 


99.99 


Mb (disk) 


Lhc<+[NII] 


90.1 


95.8 


95.8 


Mb (total) 


Lhc«+[NII] 


99.92 


99.91 


99.90 


Lhc«+[NII] 


P2cm 




99.55 


99.90 


L[oi] 


P2cm 




99.99 


99.99 


Mmdo 


P2cm 




99.99 


99.99 


Ellipticals only Mb (bulge) 


LHa+[NII] 


94.7 


79.0 


70.0 


SO and later only Mb (bulge) 


LHa+[NII] 


99.8 


99.99 


99.99 


LINERs only AVJ„t 


FWHM([NII]) 


99.69 


99.89 


99.90 


Radio-detected LINERs only AVrot 


FWHM([Nn]) 


81.5 


95.7 


95.8 


Transition nuclei only AVrot 


FWHM([Nn]) 


97.0 


98.6 


98.5 



Note. — Columns are: (1) the nuclei considered; (2) the independent variable; (3) the dependent variable; 
(4) result of the Cox proportional hazard test; (5) result of the generalized Kendall t test; (6) result of the 
Spearman p test; the values listed in columns (4), (5) and (6) are the percentage 'significance' of the correlation. 
The probability of getting the correlation inferred from the data when no actual correlation exists between the 
two variables is 100% minus the quoted value of the significance i.e. high percentages imply a low chance 
that the observed correlation is spurious. For easier reading, probabilities less than 99% arc listed to only one 
decimal place. Asurv statistical tests take into account upper or lower limits in the data. The Cox proportional 
test allows upper Umits in only one of the variables; this test could not be applied in some of the cases above. 
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Table 5 

RESULTS OF BIVARIATE STATISTICAL 
COMPARISONS FOR THE MDO SAMPLE 



Variable Variable 


Cox 


Kendall 


Spearman 


A B 




r 


p 


(1) (2) 


(3) 


(4) 


(5) 


All (154) galaxies 








MmDO P2cm 


99.99 


99.99 


99.99 


LB(bulge) P2cm 


99.99 


yy.yy 


00 00 


Ellipticals only (37 galaxies) 






MmDO P2cm 


98.5 


99.48 


99.81 


LB(bulge) P2cm 


98.4 


98.5 


99.64 


so and later only (117 galaxies) 






MmDO P2cm 


99.98 


99.99 


99.99 


Lb (bulge) P2cm 


99.78 


99.99 


99.98 


Radio detections only (66 galaxies) 




MmDO P2cm 


99.99 


99.99 


99.99 


LB(bulge) P2cm 


99.99 


99.99 


99.99 



Note. — Columns 1 to 5 are equivalent to columns 
2 to 6 of Table 4 and are explained in the footnote to 
Table 4. 



